Tae 


fet sei RA os ie: 





THE JOURNAL 





OF 


CHEMICAL PHYSICS 





VoiumeE 12, Numser 12 


DecemsBer, 1944 





Bromination of Hydrocarbons. I. Photochemical and Thermal Bromination of Methane 
and Methyl Bromine. Carbon-Hydrogen Bond Strength in Methane 


G. B. KistiAkowsky* AND E. R. VAN ARTSDALEN** 
Department of Chemistry, Harvard University, Cambridge, Massachusetts 


‘(Received October 12, 1944) 


The photochemical bromination of methane was studied in the temperature range 423-503°K 
and found to proceed through the following chain mechanism: 


(1) Bre+hvy=Br+Br 

(2) Br+CH,=CH:+HBr 
(3) CH3+Bre=CH;Br+Br 
(4) CH;+HBr=CH,+Br 
(5) Br+Br+M=Br.4+ M. 


Bromination of methyl bromide is analogous and from 7.5 to 10 times more rapid in this 
temperature range. Hydrogen bromide inhibits bromination of methane but not of methyl 
bromide. Thermal bromination was studied at 570°K and found to follow the same mechanism 
as photochemical reaction, except that bromine atoms are produced thermally. The activation 
energy of photochemical bromination of methane is 17.8 kcal./mole and that of methyl bromide 
is 15.6 kcal./mole. Varying efficiencies of different molecules as third bodies in the homo- 
geneous recombination of bromine atoms are discussed. Configurations of activated complexes 
have been assigned and by statistical mechanical calculations shown to be reasonable. Activa- 
tion energies and other data have been combined to arrive at a value for the C—H bound 
strength in methane of 102 kcal./mole at room temperature. 





UBSTITUTION reactions of bromine with 
paraffins are known to yield bromo-paraffins 
and hydrogen bromide. They have not been 
studied kinetically heretofore. Coehn and Cordes! 
have investigated the photochemical chlorination 
of methane and reported that it is analogous to 
the hydrogen chloride reaction. Other writers? 
have suggested that bromination should take 
* On leave. ; 
** Department of Chemistry, Lafayette College, Easton, 
Pennsylvania: at present on leave of absence to Depart- 
ment of Physiological Chemistry, Johns Hopkins School 
of Medicine, Baltimore, Maryland. 
(1930) and H. Cordes, Zeits. f. physik. Chemie B9, 1 


*H. J. Schumacher, X. Congr. Int. di Chimica IV, 464 
(1938); J. W. T. Spinks, Chem. Rev. 26, 129 (1940). 
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place in a similar type of chain mechanism. It was 
felt that a thorough study of the bromination of 
paraffins might lead to valuable methods for 
determining strengths of several carbon-hydro- 
gen and carbon-carbon bonds. This has proved 
to be the case.’ 


REAGENTS AND MATERIALS 


Methane of extremely high purity was pre- 
pared jointly with Dr. E. E. Roper by decom- 
position of a methyl grignard reagent with 

3 Andersen, Kistiakowsky, and Van Artsdalen, J. Chem. 
Phys. 10, 305 (1942). E. R. Van Aritsdalen, ibid. 10, 653 


(1942); H. C. Andersen and G. B. Kistiakowsky, ibid. 11, 
6 (1943). 
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ethanol. It was anticipated that oxygen would 
inhibit photochemical bromination, hence great 
pains were taken to exclude air, preparation being 
accomplished under vacuum conditions using 
reagents repurified especially. The methane was 
distilled from a trap in liquid air to one in liquid 
nitrogen (—188° to —197°) and showed a dis- 
tillation pressure between 12 and 12.5 mm. 
Vapor pressure determinations with different 
amounts vaporized indicated that the methane 
was extremely pure. 

Mallinckrodt bromine, A.R. grade, was purified 
further by storing over pure potassium bromide 
to eliminate chlorine. It then was distilled under 
vacuum through a drying tube containing P2Os. 
Alternate freezing and melting of the sample 
under high vacuum was carried out to eliminate 
traces of air. In a final distillation the middle 
third was retained for this investigation. 

A sample of methyl bromide, about 99.8 
percent pure, was donated for this work by The 
Dow Chemical Company. It was purified further 
by repeated bulb-to-bulb distillations under 
vacuum until no further change in its vapor 
pressure in a dry-ice acetone bath could be 
detected. Its vapor pressure at 0°C was 677 mm. 
Approximately 95 percent was distilled from a 
dry-ice bulb into a reservoir flask. The vapor 
pressure of the 5 percent residuum was 674 mm. 
Kemp and Egan‘ give 675.8 mm as the vapor 
pressure at 0°C. 

Hydrogen bromide was prepared by the direct 
union of electrolytic hydrogen and A.R. grade 
bromine, according to the method outlined by 
Ruhoff, Burnett, and Reid.’ Repeated bulb-to- 
bulb distillations (dry-ice to liquid air) were 
made of the product under vacuum, voiding 
considerable hydrogen bromide. The _ purified 
material formed a clear, colorless liquid and a 
perfectly white solid. Its vapor pressure was 
320 mm at —81°C. 


EXPERIMENTAL METHOD 


Bromination was followed by observing the 
rate of disappearance of bromine from the re- 
action system. For this purpose a photoelectric 
photometer of balanced construction was used. 

4J. D. Kemp and C. J. Egan, J. Am. Chem. Soc. 60, 2097 
(1938). 


5 Ruhoff, Burnett, and Reid, Organic Syntheses (John 
Wiley & Sons, Inc., New York, 1935), Vol. XV, p. 35. 
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Light of wave-length greater than 4360A was 
climinated by Corning violet filter, No. 511. 
Immediately prior to each set of measurements, 
the photometer was calibrated against known 
pressures of bromine measured with a quartz 
spiral null-point manometer against a glass scale 
mercury manometer. In general it was possible 
to reproduce bromine partial pressures to +0.1 
mm. Ability to rotate the photometer between 
two fixed paths, one of which passed through the 
reaction cell and the other through a comparison 
cell, made it possible to check the zero setting 
of the photometer immediately before and after 
each individual determination of bromine con- 
centration. 

A Pyrex glass reaction cell of about 82-cc 
volume and 106-cm? wall area was supported in 
a fixed position in an electrically heated glass 
air bath which was made large in comparison 
with the cell in order to minimize thermal 
gradients. Temperature was measured to 0.1° 
with a thermocouple and controlled to +0.2°- 
0.3°. Dead space was about 2-2.5 cc which is not 
serious for these bromination reactions in which 
there is no volume change. 

The reaction mixture was illuminated in the 
photo experiments by blue light (from a 750-watt 
Mazda projection lamp) filtered through 0.5” of 
water and 4.54-mm Corning blue-green filter 
No. 430 to eliminate red and infra-red radiation. 
Constant amperage through the lamp circuit was 
maintained by manual rheostat control. The 
axis of illumination was perpendicular to the 
photometer path. All lamps, air bath windows, 
and filters were cleaned with alcohol prior to 
each set of measurements. 

Bromine predissociation in the band region 
yields two normal bromine atoms, while dis- 
sociation in the continuum gives one normal 
atom and one excited atom. Jost® has shown 
that there is no very appreciable difference 
between using light in the band region or in the 
continuum to activate photochemical formation 
of hydrogen bromide, indicating that the elec- 
tronic excitation energy is not useful in the 
reaction. Most of the radiation in our experi- 
ments was in the continuum. Bromine was the 
only substance present which absorbed any light 
of the wave-lengths used. 

®W. Jost, Zeits. f. physik. Chemie B3, 95 (1929). 
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It is known that at the temperatures employed 
in this research no thermal decomposition of 
either methane or methyl bromide occurs.? An 
a posteriori calculation indicates that dark re- 
action (thermal bromination) should be undetect- 
able in the course of an hour at 483°K. Experi- 
mentally no dark reaction at this temperature 
was observed over a period of half an hour, in 
which time photochemical bromination was prac- 
tically complete. Therefore the photoreaction 
was stopped by inserting a metal shield between 
light source and reaction cell while bromine 
concentrations were determined with the pho- 
tometer. In this way repeated measurements of 
bromine concentration could be made after each 
desired interval of illumination. In the case of 
thermal brominations photometer readings had 
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Fic. 1. Photochemical bromination of methane. 483°K. 
Initial pressures, CH, 641 mm; Brz 45.0 mm. 


to be taken “‘on the fly,” that is, while the reading 
was changing. In all cases reproducibility was 
satisfactory. 

Bromine concentrations were measured at 
suitable times and from them rates of reaction 
during small time intervals were determined and 
an extrapolation made to obtain the rate at 
zero time. Rate constants were calculated from 
this rate at zero time. This method was employed 
because at the start of the reaction there are no 
products of reaction present to complicate the 
kinetics. It is known that all of the bromo-. 
methanes can be produced by photobromination, 
hence it is essential in a kinetic study to measure 
the rate of formation of a particular compound. 
The plot of a typical experiment is given in Fig. 1. 





_,’ H. P. Meissner and H. J. Schumacher, Zeits. f. physik. 
Chemie A185, 435 (1940). ’ 
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EXPERIMENTAL RESULTS 


It was found experimentally that the rate of 
photobromination of methane was proportional to 
methane concentration, to the square root of bro- 
mine concentration, and inversely proportional 
to the square root of total pressure. In analogy 
with the hydrogen bromide synthesis it is plau- 
sible to expect the rate determining step in 
methane bromination to involve the reaction of 
a bromine atom, produced either thermally or 
photochemically, with a methane molecule. On 
this basis the initial rate should be proportional 
to methane concentration and to the square root 
of bromine concentration, or in the case of the 
photoreaction, to the square root of absorbed 
light intensity. The Lambert-Beer law is known 
to apply quite accurately in the case of bromine 
up to a partial pressure of 50 mm,* about the 
maximum employed here, hence absorbed light 
intensity is proportional to bromine concentra- 
tion. The inverse proportionality with total 
pressure in the photochemical bromination is in 
accord with a three-body mechanism for de- 
struction of bromine atoms. 

Table I gives values of initial rate constants 
(k’) for photobromination calculated by means 
of the equation 


—d{Br2]/dt=k’[CH,][Bre}[1/P]}. (A) 


k’ is a rate constant multiplied by absorption 
coefficient and an apparatus constant. P is the 
total pressure. The three series of experiments 
are not comparable in terms of absolute rate 
because design and adjustment of the optical 
system were altered slightly between them. 
Series A and B show clearly that the initial rate 
constant remains unchanged over a threefold 
variation in both CH, and Brz partial pressures 
as well as total pressure. However, this equation 
cannot be used to evaluate rate constants during 
any appreciable part of the reaction because of 
secondary reactions and inhibition by the hydro- 
gen bromide produced. 

A direct test of the dependence of rate on the 
intensity of incident radiation was made by 
inserting between the light source and reaction 
cell a screen calibrated to transmit 44.5 percent 


8M. Bodenstein and H. Liitkemeyer, Zeits. f. physik. 
Chemie 114, 208 (1924). 
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TABLE I. Photobromination of methane. 








Initial pressures 





(mm at T°K) 
T°K CH, Bre k’(sec.~!) X104 
Series A 503 547 53.5 7.0 
618 41.7 7.3 
613 40.5 7.3 
555 41.1 won 
553 41.1 7.0 
456 41.5 aa 
311 43.0 7.3 
284 40.5 va 
279 40.8 7.0 
210 40.9 Ps 
473 25.1 7.8 
280 21.5 7.0 
497 16.3 73 
125 
Series B 483 641 45.0 1.38 
514 44.4 1.39 
210 44.4 1.41 
1.40 
Series C 483 492 44.4 2.94 
488 45.1 3.03 
486 46.2 2.89 
486 43.6 2.53 
483 44.1 2.84 
450 45.8 2.71 
2.82 
453 424 45.7 0.816 
419 $1.1 .816 
418 45.5 .830 
418 45.4 .764 
0.807 
423 406 45.2 0.196 
400 44,9 185 
0.190 








of incident radiation. If the rate is proportional 
to the square root of absorbed light and with 
validity of the Lambert-Beer law this screen 
should reduce the velocity of bromination to 
(0.445)!=0.667 or 66.7 percent of its value in 
absence of the screen. Experimentally the rate 
was reduced to 68.3 percent of its normal value. 
This is within experimental error of 66.7 percent 
and is considered as additional proof that the 
rate depends upon the square root of absorbed 
light intensity, which latter is proportional to 
- the bromine concentration. 

It was found that hydrogen bromide, a product 
of the reaction, inhibits bromination. The follow- 
ing equation was found to apply to photobromi- 
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nation of methane in the presence of HBr, 
d[Br2] k’(CH,)(Bre}'(1/P}} 
dt k{ HBr] 
1+—__—— 
k3[ Bre | 





(B) 


Table II shows results of experiments made to 
study the effect of inhibition by HBr. Values of 
k’ were taken from experiments run in absence 
of HBr. The experiments at 310°K were of very 
long duration (two months), and were carried 
out in sealed tubes, analysis being made chemi- 
cally instead of photometrically. They are sub- 
ject to lower accuracy than those at the two 
higher temperatures. In any case it is difficult to 
determine the ratio k4/k3 with high accuracy 
because the rate is rather insensitive to relatively 
large changes in HBr concentration. 

Bodenstein and Liitkemeyer,*’ employing an 
analogous expression, found for the hydrogen- 
bromine reaction a value of k4/k3=0.10 which is 
temperature independent. It is seen in Table I] 
that at high temperature HBr inhibits bromina- 
tion of CH, more strongly than that of He and 
that this inhibition has a definite temperature 
coefficient. 

An experiment was performed at 453°K to 
determine whether oxygen exerts an inhibitory 
action. To a reaction mixture of 371 mm of 
CH, and 45.7 mm of Bro, carefully dried and 
CO:-free air was admitted corresponding to 
1.6 mm of Oe. A calculation using k’ =0.807 


TABLE II. Photobromination of methane. 
Inhibition by HBr. 








Initial pressures 





(mm at T°) 
fe CH, HBr Br2 ka/ks 
483 363 121 46.6 0.212 
385 96 47.7 .204 
423 254 138 41.6 .154 
248 116 42.1 111 
310 298 143 48.4 .089 
281 135 48.4 .064 








X10-* sec.' showed that if this behaved as an 
inert gas the initial rate of disappearance of Br2 
should have been 0.60 mm/min. The experi- 
mental value was 0.4 mm/min., indicating strong 
inhibition by oxygen. However, this inhibition 
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was present at the start of the reaction only, the 
rate having increased to 0.6 mm/min. after about 
1.5 mm of Bro had been consumed. The reaction 
was normal from then on. This experiment is in 
agreement with known behavior of oxygen in 
other brominations and bromine sensitized oxi- 
dations.’ 

Series C of Table I gives data for a 60° tem- 
perature interval which have been used to calcu- 
late the activation energy of photochemical 
bromination of methane. Before such data may 
be employed with assurance it is necessary to 
answer two questions: (1) Does the absorption 
coefficient of bromine vary appreciably over the 
temperature range in question, and (2) is the 
ideal gas law sufficiently exact under the condi- 
tions of the experiment? The second question 
was tested by applying the Berthelot equation 
to Bre. Deviations between it and the ideal gas 
law in our pressure ranges were considerably less 


than 0.5 percent, hence the ideal gas law was' 


used. Bodenstein and Liitkemeyer® have shown 
that the absorption coefficient of bromine varies 
about 5 percent between room temperature and 
490°K at 5350A, decreasing with shorter wave- 
length to 1 percent variation at 5000A. Rela- 
tively little light of wave-length greater than 
5000A was used. Consequently no correction 
was made for temperature dependence of the 
absorption coefficient over our smaller tempera- 
ture interval. In Fig. 2, log k’ has been plotted 
against 1/T giving a good straight line. The 
slope of this line gives 18.25 kcal. for the un- 
corrected energy of activation. Taking the mean 
temperature of 453°K, RT/2=0.45 kcal. leaving 
17.8 kcal. by subtraction from 18.25 for the 
corrected energy of activation.® 

The photochemical bromination of methyl 
bromide corresponds to the second step in the 
over-all bromination of methane. It was studied 
kinetically in the same manner as was methane. 
The rate of disappearance of bromine follows the 
kinetic expression 


—d[Br2]/dt=k"[CH3Br][Br2}[1/P}!. (C) 


Rate constants are given for three series of 
experiments in Table III. Series C has been 
used to evaluate the activation energy as in the 


* This correction is based on the simple bimolecular rate 
equation, R=Z'Tte-#/RT. 
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case of methane, the data being plotted in Fig. 2. 
The corrected value for the activation energy is 
15.6 kcal./mole. The effect of hydrogen bromide 
was investigated at both 483°K and 423°K. No 
inhibition was observed. 
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Fic. 2. Temperature coefficients. Each circle represents 
average of all runs at that particular temperature. 


The rate constants for bromination of methane 
and methyl bromide being known separately, it is 
possible to calculate what the rate should be in 
mixtures of methane and methyl bromide. Satis- 
factory agreement was obtained between experi- 
ment and calculation when tested. 

Several experiments at both 483 and 423°K 
were made on the photochemical hydrogen- 
bromine reaction in order to check our method 
on a very well studied reaction. Good repro- 
ducibility was obtained. The activation energy 
calculated from these experiments was 17.5 
kcal./mole in agreement with Bodenstein’s value 
of 17.6 kcal./mole. The data are plotted in 
Fig. 2. 

The rates of thermal bromination of methane, 
methyl bromide and hydrogen were measured at 
570°K. The method of measurement was essen- 
tially the same as that used in the photochemical 
experiments. Data are recorded in Table IV. 
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Rate constants were calculated from the initial 
rate of reaction using the equation: 


—d[Bre ]/dt=kCCH, ][Bre }}. (D) 


It is observed that the rate is independent of 
total pressure and that concordant values of k are 
obtained over a twofold variation of partial 
pressures. 

The data of Bodenstein and Lind!® were used 
to calculate a rate constant for the hydrogen 
reaction at 570°K of k=0.022 which is in good 
agreement with our value of k=0.0196. A differ- 
ence of 1 to 1.5° in the temperature scales would 
bring the two values into excellent agreement. 


DISCUSSION 


The following chain is postulated as the re- 
action mechanism operative in bromination of 
methane. 


Bre+hv—Br+Br (1) 

or Br2+M—Br+Br+M (la) 
Br+CH,—-CH;3;+HBr (2) 
CH3;+Br.2—-CH3Br+ Br (3) 
CH3;+HBr—CH.+ Br (4) 


Br+Br+M—-Br:+M (5) 


Br+CH;Br—-CH.Br+HBr (2a) 


CH2Br+Br2—CH2Br2+ Br (3a) 
CH.Br+HBr—CH;3Br+Br (4a) 
etc. 


At sufficiently low pressures one would expect 
chain termination to result from diffusion of 
bromine atoms to the wall followed by hetero- 
geneous recombination. This is not the case in 
our experiments, for the rate then would be 
first order with respect to Bre concentration 
instead of 0.5 order. 

Reactions involving production of hydrogen 
atoms have been ruled out on the basis of energetic 
considerations. Some would involve Walden type 
inversions which are known to have high activa- 
tion energies. For example, it has been shown" 
that the reaction CH,+H—CH;+H2: requires 8 

10 M. Bodenstein and S. C. Lind, Zeits. f. physik. Chemie 
57, 168 (1906). 


1 Gorin, Kauzmann, Walter, and Eyring, J. Chem. Phys. 
7, 633 (1939). 
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TABLE III. Photobromination of methyl bromide. 


Initial pressures 
(mm at T°K) 








T°K CHs3Br Br2 k’’(sec.~!) K 108 

Series A 503 337 52.4 5.36 
159 41.8 5.15 
347 40.6 5.19 
Series B 483 381 43.2 1.07 
453 42.8 1.02 
398 42.6 1.08 
298 32.9 Lae 
211 29.0 1.17 
Series C 483 259 48.5 2.07 
287 45.9 2.21 
255 45.2 2.25 
472 44.1 1.93 
2.12 

453 331 45.5 0.754 

336 44.5 .769 

0.762 

423 317 45.0 0.209 

308 44.4 .199 

295 42.8 .198 

0.202 





TABLE IV. Thermal brominations at 570°K. 





Initial pressures 
(mm at 570°K) 
X 








Bre k(ce mole~!)3 sec.~! 

A. —CH, 517 53.0 0.0247 
331 55.5 .0250 

317 48.1 .0249 

227 51.0 .0240 

353 25.0 .0253 

325 22.0 .0247 

0.0248 

B. —He 431 51.2 0.0193 
405 52.0 .0198 

0.0196 

C. —CH;Br 352 52.7 0.156 
210 $3.2 .163 

342 25.6 .162 

0.160 


to 11 keal. activation, energy while the Walden 
inversion CH,+H—-CH.i+H requires 35 to 40 
kcal. Thus there is good evidence against re- 
actions of the type CH4+ Br—-CH;Br+H. 
Using the above reaction mechanism and 
assuming steady state concentrations of chain- 
propagating radicals, one can derive the ex- 
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pression 
d{CHsBr] ke[CHs](n/ksLM)) 
dt ; k(HBr] 7 
k3[ Bre ] 


where 1 represents the number of absorbed light 
quanta and [17] is the total concentration of all 
molecules present. It is clear that 1=k,~,al( Bre |, 
where ka» is a temperature independent appa- 
ratus constant, J is the intensity of incident 
radiation, and a@ is the absorption coefficient of 
bromine. Because J was a constant in any set of 
experiments we may write k’=ko(Rapal/ks)?. 
Hence, Eq. (E) is seen to be equivalent to 
Eq. (B). It is based on methyl bromide as the 
only bromomethane formed and is, therefore, 
a first-order approximation only. Considering 
secondary bromination to produce methylene 
bromide with negligible inhibition of this step 


TABLE V. Frequency factor and entropy of 
activation at 570°K. 





A x S 
Reaction Expt. Collision theory Expt. Calc. 
Br+He 1.14 K10" 11.7 X10" —2.31 
Br+CHy, 1.48 X10" 7.2 X10!2 —1.81 —1.84 
Br+CHs3Br 1.35 X10" 5.9 K10" —1.98 —1.94 


by hydrogen bromide, the second-order approxi- 
mation for the rate becomes 


d(CH;Br]  k’([CH,][Br2}'[1/P]} 
4 AHA 
k3[ Bre | 
—k"(CH;Br][Bro}{{1/P]}, (F) 





where k’’ = koa(kapal/ks)'. The rate of disappear- 
ance of bromine clearly will be the sum of the 
two terms on the right. This equation was found 
to fit our experimental data when applied by 
means of successive approximations. 

From the postulated mechanism one derives 
the first-order approximation for the thermal 
reaction : , 


d(CH;Br ] kLCH, j[Brz }} 
dt - k[HBr] ’ 


k3[ Bre | 





(G) 
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in which k=ke(Keq)! where Keq is the equi- 
librium constant for dissociation of bromine. 
This reduces to Eq. (D) at zero time when the 
concentration of hydrogen bromide is zero. 

Qualitatively we may observe that the rate of 
disappearance of bromine will increase for a time 
after the reaction begins because of the gradual 
increase in concentration of methyl bromide 
which is brominated from 7.5 to 10 times more 
rapidly than methane in the temperature range 
we investigated. Inhibition by hydrogen bromide 
becomes progressively more pronounced in the 
latter stages of bromination when its concentra- 
tion becomes relatively high in respect to bro- 
mine. Acceleration by methyl bromide is more 
important than inhibition by hydrogen bromide. 
However, decreasing bromine concentration of 
itself also causes a decreasing rate of bromina- 
tion. These effects are seen in Fig. 1. 

It is evident from the form of the kinetic 
equation that the rate determining step is re- 
action (2) in the case of methane and (2a) in 
the case of methyl bromide. If this were not 
true the rate would depend on some other con- 
centration functions. Energetic considerations 
also are in agreement. Reaction (3) will be 
exothermic to an appreciable extent because the 
C—Br bond is stronger than the Br—Br bond 
(45.2 kcal.) and because it is a free radical re- 
action it should require very low or zero activa- 
tion energy. On the other hand, reaction (2) 
should be endothermic because the C—H bond 
is certainly stronger than the H—Br bond (85.6 
kcal.). It must, therefore, require an activation 
energy at least equal to its endothermicity. On 
the basis of these considerations the activation 
energy of bromination of methane, 17.8 kcal., 
has been ascribed to reaction (2), while the 
corresponding value of methyl bromide, 15.6 
kcal., has been assigned to reaction (2a). 

The temperature dependence of inhibition 
by hydrogen bromide enables a determination 
of the difference in activation energies between 
reactions (3) and (4) and a reliable estimate 
of the activation energy of reaction (4). Em- 
ploying the usual, simple bimolecular equation 
k=AT'e-=/®T, where A is a temperature inde- 
pendent frequency factor and E the activation 
energy, we express the logarithmic ratio of the 
rate constants of reactions (3) and (4) in the 
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TABLE VI. Rate constant ratios. 
kcu/Fu: RCH;Br/*H kCHsBr/FCH, 
Calc'd with Calc’d with Calc'd with 
Experi- efficiency Experi- efficiency Experi- efficiency 
T°K mental i314 34:22 mental $33 260222 mental 21 aa 
Thermal 1.26 8.2 6.4 
570 
Photo 
483 1.39 1.05 1.35 10.5 9.6 10.5 iS 9.0 7 
453 9.3 10.9 9.4 
423 1.33 1.01 1.30 14.2 12.8 14.0 10.6 [2:5 10.7 
form H2=2.76A, CH4=3.2A, CH3Br=4.5A, Br=3.0A. 
A 4. EimE Values of A and Z’ (= Z/T*) are given in Table V. 
Sea FS pectin It is interesting that A and Z’ differ little from 
ks As RT one another, the steric factor being approxi- 


For these two similar reactions the assumption 
is made that the quantity (E,—E3) does not 
vary appreciably with temperature. Since A also 
is essentially temperature independent, the ex- 
pression can be differentiated to give 


ky Es—E; 
dia —o——_—41 


ks RT? 


which combined with inhibition data (Table 11) 
yields (E,—E;)=2 kcal. For reasons already 
mentioned it is not unreasonable to assume that 
the activation energy of reaction 3 is zero, 
making 2 kcal. the activation energy of reaction 4 
in good agreement with the value of 1.5 kcal. 
obtained by Andersen and Kistiakowsky” from 
experiments on the photolysis of methyl iodide 
in the presence of hydrogen bromide. 

Absolute rates of thermal bromination at 
570°K and data of Gordon and Barnes on the 
thermal dissociation of molecular bromine have 
been used to calculate the temperature inde- 
pendent frequency factor A in the equation 
k=AT'e-=/8T as applied to the rate determining 
single reactions of bromine atoms with hydrogen, 
methane, and methyl bromide. Experimental 
values of E determined in the photochemical 
brominations were used. The collision number 
Z was calculated by classical collision theory 
employing the following molecular diameters, 


2H. C. Andersen and G. B. Kistiakowsky, J. Chem. 
Phys. 11, 6 (1943). 

13 A. P. Gordon and C. Barnes, J. Chem. Phys. 1, 692 
(1933). 


mately 0.1 to 0.2. This is in contrast to the 
reaction of hydrogen atoms with methane where 
a steric factor of the order of 10- to 107% is 
required." A possible explanation of this differ- 
ence may be found by considering the relative 
sizes of bromine and hydrogen atoms. The 
hydrogen atom presumably must approach the 
methane molecule more precisely oriented along 
the axis of a C—H bond than is required in the 
case of the bromine atom. 

Experimental entropies of activation were cal- 
culated from the data of thermal experiments at 
570°K using a formula of Kistiakowsky and 
Ransom," 


kpT 


k=—e’eS/Re—E/RT 


oh 


where S is the entropy of activation, kg the 
Boltzmann constant, and h Planck’s constant. 
“Theoretical” entropies of activation were calcu- 
lated by statistical mechanical methods. Con- 
figurations of activated complexes were assigned 
by analogy with those calculated by Gorin, 
Kauzmann, Walter, and Eyring" for reaction 
between hydrogen atoms and methane. It was 
assumed that the minimum activation energy is 
required when a bromine atom approaches a 
methane molecule along the axis of a C—H 
bond. No change in configuration of the CH; 
part of the molecule was assumed. Interatomic 
distances and vibrational frequencies were as- 


4G. B. Kistiakowsky and W. W. Ransom, J. Chem. 


Phys. 7, 725 (1939). 
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signed as follows, 


H;—-—_C------ H------------Br 
109A 1.17A 1.86A 
C—H stretching 3000 cm=! (4) 
CH; torsional 1440 (6) 
H—Br stretching 2200 (1) 
Br—H—CH,; bending 500 (2). 


In the case of methyl bromide the assignments 
were 


Bri.91A 117A 1.86 


Frequencies of methyl bromide’ plus others for 
bromine atom assignments were 3061 (2), 2900 
(1), 2200 (1), 1450 (2), 1305 (1), 957 (2), 610 (2), 
and 105 (2). A comparison of experimental and 
“theoretical” entropies of activation is given in 
Table V. It is observed that these complexes 
have low, doubly degenerate bending frequencies 
which make large contributions to the entropy. 
Unfortunately they cannot be estimated with 
high precision. However, agreement between 
experimental and calculated entropies of activa- 
tion indicates that the structures assigned to 
the activated complexes are reasonable. 

Ratios of the experimental rate constants for 
thermal and photochemical reactions are shown 
in Table VI. Clearly in. the case of thermal 
bromination these are ratios of the rate con- 
stants of reaction (2) and its prototypes, the 
equilibrium constant for dissociation of bromine 
cancelling out. However, for photochemical bro- 
mination they would represent the ratios of the 
bromine atom reactions only if all molecules were 
equally effective third bodies for homogeneous 
recombination of bromine atoms. That such is 
not the case is shown in Table VI. Ratios for 
photobromination calculated from thermal abso- 
lute rates and experimental activation energies 
on the assumption of equal efficiencies do not 





* D. P. Stevenson and J. Y. Beach, J. Chem. Phys. 6, 
25 (1938). 
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agree with experiment. Rabinowitch and Wood? 
have studied the influence of various molecules 
on recombination of bromine atoms and report 
that the efficiencies of hydrogen and methane 
are in the ratio 2.2:3.6, respectively. When these 
values are used in calculating rate constant 
ratios excellent agreement with experiment is 
obtained as shown in Table VI. We deduced from 
our experiments that the efficiency of methyl 
bromide on the same basis as hydrogen and 
methane was 2.7, the results being indicated in 
the table. It is surprising that methyl bromide 
efficiency should be less than that of methane. 
However, this quantity is very susceptible to a 
small difference in the values of activation 
energies. A recalculation assuming an activation 
energy of 16.2 kcal., i.e., 0.6 kcal. above our 
accepted value of 15.6 kcal., for bromination of 
methyl bromide yields an efficiency of 3.9. It can 
be said that our experiments confirm the work of 
Rabinowitch and Wood and indicate that the 
efficiencies of methyl bromide and methane are 
roughly equal. In the inhibition experiments with 
hydrogen bromide no account has been taken of 
varying efficiency of methane and hydrogen 
bromide. Work of Hilferding and Steiner!’ com- 
bined with that of Rabinowitch and Wood indi- 
cates that the efficiencies of methane and hydro- 
gen bromide are equal. Differing efficiencies 
would change the absolute value of the inhibition 
constant but should not affect its temperature 
coefficient. 

The experimental rate constant ratios of photo- 
chemical bromination in Table VI can be used as 
an estimate of the relative quantum yields. If the 
quantum yield for hydrogen bromide synthesis 
is one'® then it is about 1.4 for formation of 
methyl bromide. 


CARBON-HYDROGEN BOND STRENGTH 


Values for activation energies of reactions (2) 
and (4) may be combined to obtain the heat of 
reaction at the mean experimental temperature 


16E. Rabinowitch and W. C. Wood, Trans. Faraday 
Soc. 32, 907 (1936). 

17K. Hilferding and W. Steiner, Zeits. f. physik. Chemie 
B30, 399 (1935). 

18 See, for example, W. A. Noyes, Jr., and P. A. Leighton, 
The Photochemistry of Gases (Reinhold Publishing Corpo- 
ration, New York, 1941). 
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453°K for CH4+Br—CH;+HBr. Because the 
numbers of molecules of products and reactants 
are equal, we may write 


AH 453 = AE 453 = (Ex— Ex) = (17.8— 2) = 15.8 kcal. 
AH,? was calculated using the following values of 


0 0. 
at. 
ott, 


CH, 3.83 kcal. /mole 
CH; 3.75 kcal./mole 
HBr 3.15 kcal./mole 
Br 2.25 kcal./mole 
H 2.25 kcal./mole. 


Values for CH, and HBr are from a review by 
Pitzer,’® H and Br were calculated for free atoms 
with translational degree of freedom only, and the 
value for CH; is that for CH, diminished by 0.3 
times the vibrational contribution. We thus ob- 
tain AH,»°=15.0 kcal. This may be combined 
with the value 85.8 kcal. for the heat of dissocia- 
tion of hydrogen bromide” to obtain the carbon- 
hydrogen bond strength as follows: 


CH,+Br—-CH;+HBr 
HBr—-H+Br 


AH,?= 15.0 kcal. 
AH,°= 85.8 kcal. 


AH,®= 100.8 kcal. 





CH.—-CH;+H 


The value of A/7® at any temperature may be 
calculated in similar manner. We obtain the 
following rounded values for the carbon-hydrogen 





i9 K, Pitzer, Chem. Rev. 27, 39 (1940). 

20 D(H—Br) =85.8 kcal. by combination of the heat of 
dissociation of He = 102.48 kcal. [H. Beutler, Zeits. f. physik. 
Chemie B29, 315 (1935) ] with that for Brz=45.22 kcal. 
and heat of formation of HBr from its elements = 11.99 
kcal. (reference 13). 
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bond strength at the indicated temperatures. 


T°K D(CH3;—H) 
0 101 kcal. 

293 102 kcal. 

453 103 kcal. 


Although quite divergent values of the carbon- 
hydrogen bond strength in methane have been 
reported, a considerable body of evidence points 
to its being in the range 100 to 104 kcal./mole.”! 
The agreement between our value and that of 
Stevenson”'£ derived from electron impact critical 
potentials is particularly impressive. Thus two 
entirely different experimental methods of deter- 
mining the carbon-hydrogen bond strength agree 
within their experimental errors (about 1 kcal.). 
Semi-empirical calculations of H. A. Taylor and 
Burton*!* by the Voge method yield the value of 
102.1 kcal. if one assumes the zero-point energy 
of CH; to be three-quarters that of methane. 

Similar assumptions regarding activation ener- 
gies lead to the conclusion that the substitution 
of bromine for one hydrogen in methane lowers 
the bond strength by about 2 to 3 kcal. von 
Braunwarth and Schumacher” found an activa- 
tion energy of the order of 10 kcal. for bro- 
mination of bromoform, which indicates that the 
introduction of 3 bromine atoms into methane 
lowers its C—H bond strength by about 7-8 kcal. 
Probably the introduction of two bromine atoms 
lowers this bond strength about 5-6 kcal. 


*t (a) L. Pauling and J. Sherman, J. Chem. Phys. 1, 606 
(1933), 100 keal.; (b) R. G. W. Norrish, Trans. Faraday 
Soc. 30, 103 (1934), 103 kcal.; (c) F. O. Rice and K. K. 
Rice, Aliphatic Free Radicals (Johns Hopkins Press, Balti- 
more, 1935), p. 75, 100 kcal.; (d) F. Patat, Zeits. f. physik. 
Chemie B32, 290 (1936), 102.5 kcal.; (e) H. A. Taylor and 
M. Burton, J. Chem. Phys. 7, 572 (1939), 102.1-103.7; (f) 
E. C. Baughan and M. Polanyi, Nature 146, 685 (1940), 
103.6 kcal.; (g) D. P. Stevenson, J. Chem. Phys. 10, 291 
(1942)—101 keal. 

22 yon Braunwarth and H. J. Schumacher, Kolloid Zeits. 
89, 184 (1939). 
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Ethyl Bromide. Carbon-Hydrogen Bond Strength in Ethane 
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The kinetics of the photochemical bromination of ethane was determined. The mechanism 
is very similar to the analogous brominations of hydrogen and methane: 
(1) Bro+hy=Br+Br 
(2) Br+C:He=C»H;+HBr 
(3) CsH;+Bre=C2.H;Br+Br 
(4) CoH; +HBr=C.H,+Br 
(5) Br+Br+M=Bro4+ M. 


Rate constants computed ‘by assuming various efficiencies for ethane as the third body in 
step (5) are presented, as well as those calculated by considering heterogeneous recombination 
of bromine atoms. The rate of bromination of ethyl bromide is about the same as that for 
ethane, but its exact determination was precluded by the presence of rapidly brominated 
impurities. The experimental activation energy, 13.6+0.5 kcal., is quite insensitive to the 
several variations in mechanism discussed, and indicates a value of 98+2 kcal. for the upper 
limit of the strength of the C2H;—H bond. This figure is in good agreement with absolute values 


reported by other workers. 


INETIC studies of the bromination of meth- 
ane! having proved successful, it was de- 
cided to extend the method to ethane. 


EXPERIMENTAL 


thane was prepared by catalytic hydrogena- 
tion with cobalt-nickel catalyst of 99.5 percent 
ethylene, a standard product of the Ohio Chemi- 
cal Company. The hydrogen was purified care- 
fully to remove traces of oxygen and passed 
through in twofold excess. Any unsaturates re- 
maining after the hydrogenation were removed 
by fuming sulfuric acid; potassium hydroxide 
was used to free the ethane from acid gases. 
The ethane was condensed in liquid nitrogen. It 
was then given a batch distillation with efficient 
stirring followed by several bulb-to-bulb dis- 
tillations under high vacuum to remove perma- 
nent gases and passed through a glass wool filter 
at dry ice temperature to remove mercury. 
Purity tests described by Nickle,? indicated about 
. i ene address: Kellex Corporation, New York, New 

ts Department of Chemistry, Lafayette College, Easton, 
Pennsylvania. At present on leave of absence to Depart- 
ment of Physiological Chemistry, Johns Hopkins School 
of Medicine, Baltimore, Maryland. 

'G. B. Kistiakowsky and E. R. Van Artsdalen, J. Chem. 


Phys. 12, 469 (1944), hereafter termed I. 
* A. G. Nickle, Thesis (Harvard University, 1941). 
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0.01 percent unsaturates, as ethylene, and 0.02 
mole percent high boiling material, calculated as 
butane. The ethane registered an equilibrium 
pressure of 0.02 mm Hg when immersed in liquid 
air. (Compare Kistiakowsky and Nazmi: 0.95 
xX 10-* at liquid oxygen temperature.) 

Two samples of ethyl bromide were used. The 
first sample was prepared from Baker’s ‘‘purified”’ 
ethyl bromide by distilling through a 9-foot 
glass-helix packed column and recovering the 
fraction boiling at 38.6-38.8°C. This fraction 
was further purified by three bulb-to-bulb high 
vacuum distillations. It had to be stored behind 
a stopcock lubricated with Shawinigan low melt- 
ing resin because of solubility in ordinary grease. 
The second sample was prepared from the same 
“Purified” grade by treatment with liquid bro- 
mine for one-half hour, subsequent removal of 
bromine by 20 percent sodium hydroxide and 
drying of the product over calcium bromide. 
The dried material was distilled as above; a frac- 
tion boiling at 38.4—38.5° being redistilled and a 
middle cut boiling at 38.3° (758 mm corr.) taken 
for final bulb-to-bulb vacuum distillation. It was 
stored behind a glass valve. It will be shown that 


3G. B. Kistiakowsky and F. Nazmi, J. Chem. Phys. 6, 18 
(1938). 









































480 H. C. ANDERSEN AND 
these samples of ethyl bromide apparently con- 
tained appreciable amounts of impurity, prob- 
ably ethanol. 

Bromine and hydrogen bromide were prepared 
and purified as described in I. 

The experimental method was essentially the 
same as described in 1. However, one important 
change was made. Small pressures were measured 
with the quartz spiral null-point manometer 
against a butyl-phthalate manometer instead of 
a mercury manometer. Special pains were taken 
to outgas the butyl-phthalate thoroughly before 
each set of measurements. Butyl-phthalate pres- 
sures were converted to mm of mercury employ- 
ing measured absolute densities of the butyl- 
phthalate over the temperature range 20-31°C. 

In the kinetic experiments, partial pressures 
of bromine were measured after various periods 
of illumination. The rate of reaction was deter- 
mined over the first ten percent of reaction only, 
during which interval interference by products of 
reaction is at a minimum. It was found that the 
rate during the first ten percent was practically 
linear. Constants in this paper properly are 
termed initial constants in the same way as 
those cited in I. 


RESULTS 


Bromination of ethane conceivably could lead 
to a considerably greater diversity of products 
than in the case of bromination of methane. 
It was necessary to prove the assumption that 
during the first 10 percent of reaction essentially 
only ethyl bromide is formed by photo-bromina- 
tion of ethane. For this purpose a 2-liter mixture 
of 43.5 mm of bromine and 495 mm of ethane 
was illuminated until reaction was practically 
complete and an analysis made of the reaction 
products. From average determinations of vapor 
pressure and refractive index it was concluded 
that the product of reaction consisted of a 
minimum of 65 percent of ethyl bromide. Under 
the experimental conditions employed this corre- 
sponds to a minimum of 97 percent of ethyl 
bromide during the first 10 percent of reaction. 
It was on this basis that constants were calcu- 
lated from rates during the initial 10 percent of 
reaction. 

The measured thermal, or ‘‘dark,’’ rate of 











BE. KR. VAN ARTSDALEN 
TABLE I. Photobromination of ethane. 
Initial pressure Measured 
Temp (mm at (°C) rate k k’ 
(°C) Bre CoHe (mm/min.) (sec.~! K 105) 
(corrected to 89.8°) 
90.1 17.40 660.0 10.4 161 44 174 
89.6 17.38 608.8 9.3 155 43 172 
90.2 17.64 319.6 5.92 132 42 174 
89.4 16.52 317.1 6.11 147 46 183 
89.9 16.35 310.0 5.69 136 43 166 
(corrected to 59.7°) 
59.9 32.7 664.9 2.18 25 7.6 26 
59.4 14.7 648.0 1.60 28 7.4 30 
59.3 15.25 733.5 1.63 26 6.9 28 
59.5 9.48 721.8 1.70 34 8.6 37 
59.7 9.52 741.4 1.99 39 9.8 42 
59.5 16.82 288.5 1.09 27 8.6 33 
59.1 16.05 288.2 .96 25 7.8 31 
59.1 15.92 287.8 1.06 28 8.7 35 
60.1 15.9 290.4 1.10 27 8.4 36 
60.2 10.82 278.5 98 29 8.5 39 
59.7 31.42 174.1 81 20 8.7 26 
58.6 29.55 188.2 sae 20 8.3 26 
59.5 16.2 187.4 .76 24 8.5 36 
59.9 9.80 192.8 72 28 8.5 47 
59.3 30.52 53.5 .226 12 7.4 31 
(corrected to 35.6°) 
35.1 14.1 272.3 17 4.7 1.5 6.3 
35.2 14.80 282.1 .221 5.9 1.88 7.4 
36.7 14.68 281.5 .231 5.6 1.77 7.5 
35.4 14.50 274.6 .213 5.8 1.82 7.8 





bromination at 96°C was 0.02 mm/min. Such 
“dark”’ reaction was negligible, since rates of 
photobromination were of the order of 5-10 
mm/min. at 90°C. 

If bromination of ethane, in its early stages, 
involves a mechanism analogous to that for 
bromination of methane, the rate of the photo- 
chemical reaction should be proportional to the 
square root of light intensity. An experiment 
performed to check this gave fair agreement with 
prediction (observed rate=0.62 mm per min.; 
calculated rate =0.70 mm per min.), but indicated 
that the rate probably depends upon a power 
of the light intensity between 0.5 and 1. 

Data for a series of experiments at three mean 
temperatures are presented in Table I. All rate 
constants in this paper have been corrected to 
mean temperatures with the aid of temperature 
coefficients calculated from the experimentally 
evaluated activation energy. The coefficients 
were 1.073 at 35°, 1.062 at 60°, and 1.052 at 90°. 
This procedure will be justified later when it will 
be shown that the activation energy is only 
slightly subject to the “‘constant”’ used in calcu- 
lating it. The corrections rarely exceeded 2 
percent. 

k was calculated from the following rate 
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equation in analogy with the methane-bromine 





reaction. 
—d[Brz | d{C2H;Br | k(C2He [Br }} 
dt 7 dt 7 [P}} 


where P=total pressure. This equation is valid 
for initial rates only, where reaction products do 
not interfere. 

Definite trends are observed in k, viz., de- 
creasing k with decreasing ethane pressure and 
increasing k with decreasing bromine pressure, 
indicating that the power of ethane concentration 
should be greater than one and that of bromine 
less than 0.5. Let us consider the following 
reaction mechanism in analogy with the bromina- 
tion of methane. 


Bre+hv—Br+Br (1) 
Br+C2H,-C:H;+HBr (2) 
C:H;+ Br2—-C2H;Br+Br (3) 
C.H;+HBr—-C,.H,.+ Br (4) 
Br+Br+ M—Br.4+ M. (5) 


Equation (A) can be derived from this mechanism 
for the condition that [HBr ]=0.k=ko(Rapal/ks)? 
where a=absorption coefficient, J is the incident 
light intensity, and k,, is a temperature inde- 
pendent apparatus constant. Two effects may be 
important in accounting for the observed trends 
in k, namely, different efficiencies of molecules in 
promoting reaction (5) and heterogeneous re- 
action. 

If the efficiency of the third body M in reaction 
(5) is not the same for all molecules in the re- 
acting system the rate of this reaction will not 
be given by ks{_M ][Br }* but rather by a summa- 
tion over individual effects of the various mo- 
lecular species, such as k;[uA+nB+--- ][Br . 
That different molecules have different efficiencies 
in this respect is well known,‘ and has been 
demonstrated in the case of the methane re- 
action.! However, it is quite general that more 
complicated molecules are better able to absorb 
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translational energy of colliding bromine atoms 
than are simpler molecules, methane being shown 
to be more effective than hydrogen. The deriva- 
tion of the rate equation is such that ethane 
would have to be considered less efficient than 
bromine to account for the observed effect. By 
assuming that each three-body collision between 
two bromine atoms and a bromine molecule is as 
effective as n such collisions in which ethane is 
the third body the following rate expression is 
derived when no HBr is present. 


—d[ Br | . k'[ Bre }'[CeHe | 
7 ( —) (B) 





[Bre ]+———— 
n 

where k’=k2'(Rapal/ks)'. The constants k’ in 
Table I have been calculated with »=20 which 
gave the best set of values. This would indicate 
that ethane is only about one twentieth as effec- 
tive a third body as bromine, which is a rather 
surprising result. It should be noted, however, 
that fairly large changes in n, e.g., m= 25 instead 
of 20, reflect only moderate changes in k’. 

Heterogeneous recombination of bromine atoms 
on the walls of the reaction vessel would ac- 
count for a decrease in rate constant with de- 
creasing total pressure. Rabinowitch and Wood‘ 
have demonstrated this effect, but in general 
such heterogeneous recombination is not im- 
portant at pressures above 200 mm, although 
subject to reaction vessel dimensions and in- 
tensity of illumination. It is possible to make an 
empirical correction for the wall reaction. We 
add to the above reaction scheme the following 
recombination reaction. 


wall 

Br+Br——Brz. (6) 
The rate of disappearance of Br atoms is in- 
creased by the quantity 2k Br ]/L1/], where ke 
is a constant which depends upon vessel dimen- 
sions, diffusion coefficient, etc. The rate of re- 
action derivable from this scheme for zero HBr 
concentration is 


+{4108n]) -—-| 7 





dt 


689 (1935). 


‘[M] 


* E. Rabinowitch and W. C. Wood, Trans. Faraday Soc. 32, 907 (1936); E. Rabinowitch and S. C. Lehmann, ibid. 31, 
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where k”’=ko!’(Rapal)?/ks? and shows the tem- 
perature dependence of bromination, ascribed 
to ke, and C?=k,?/4k,ks. We are left with an 
adjustable constant C which was evaluated from 
two experiments in which total pressure varied 
by a factor of about 8. The value found was 
C=11 which was used to compute k”’ in Table I. 
The effect of varying C was tested and found to 
be small for C=9 and 13. It is observed that this 
treatment gives good constancy for high and low 
total pressures, but there is a marked trend toward 
higher constants for decreasing bromine pressure. 

Perhaps a combination of the two modes of 
correction (varying third-body efficiency and 
wall recombination) would improve the situation, 
but the use of two adjustable parameters would 
be involved and scarcely justified. 

Constants calculated according to the three 
methods outlined above, and presented in Table 
I, have been used to compute activation energies. 
Considerations concerning computation of activa- 
tion energy of photobromination of methane 
outlined in I are applicable in this case also. 
The following results were obtained : 


from k E=13.4+0.7 kcal. 
from k’ E=13.3+0.4 kcal. 
from k” E=14.0+0.4 kcal. 


The average experimental value, E=13.6+0.5 
kcal., seems quite certain. It is noted that al- 
though absolute magnitudes of rate constants 
calculated by the several methods vary greatly, 
their temperature coefficients remain practically 
unchanged. The experimental energy of activa- 
tion has been corrected to the true value by 
subtracting the term 3RT which yields a final 
value of 13.3+0.5 kcal. Since in the proposed 
reaction scheme the temperature dependence may 
be ascribed tothe endothermicreaction (2) between 
Br atoms and CsHg molecules, we have assigned 
an activation energy of 13.3 kcal. to this reaction. 
Various experiments were made to measure the 
extent of inhibition by HBr, which would occur 
through reaction (4) Cs,H;+HBr—-C:H,.+Br. 
The rate equation in the presence of HBr is 


—d[Br2]_&[Br2]\[C:He][1/P} 
dt be (HBr 
ks [Br ] 





» (C) 


R. VAN ARTSDALEN 
where k= ko(Rapal/ks)*. In actual practice experi- 
ments with HBr were made in pairs bracketed 
alternately between those with no HBr present. 
Data at three temperatures are presented in 
Table II. It is evident that inhibition by hydro- 
gen bromide is temperature independent, in con- 
trast with the effect observed in bromination of 
methane.' Temperature independence of the ratio 
of these two specific reaction rate constants 
k4/ks indicates that the activation energies of 
the two reactions are essentially equal. The 
magnitude of k4/k3 is considerably larger than in 
hydrogen bromide synthesis® or methane bro- 
mination,! from which one may surmise that the 
steric factors of the two reactions are more nearly 
equal in the case of ethane than in the case of 
hydrogen or methane. 

A- few experiments at 60-61°C were made of 


TABLE II. Photobromination of ethane. Inhibition by HBr. 





Initial pressures (mm at (°C) 
“oH HE 


°c Co 3r Bre ka/ks 
90.0 326.3 102.3 16.2 0.68 
89.9 S25. 103.0 16.97 74 
89.8 323.1 103.9 16.39 68 
89.9 312.9 104.0 16.2 58 
89.9 309.5 103.5 16.6 58 
60.6 195.7 98.5 15.40 64 
35.2 266.0 82.8 14.2 7 





the photobromination of ethyl bromide. It was 
found that there was an extremely rapid initial 
reaction in which approximately 7 percent of the 
ethyl bromide appeared to react. This was 
followed by a nearly linear reaction taking place 
at a rate about equal to the rate of bromination 
of ethane under similar conditions. Upon re- 
admitting fresh ethyl bromide to the reaction 
mixture a second rapid reaction corresponding 
to about 7 percent of the readmitted ethy! 
bromide occurred. The initial rapid reaction was 
shown also by the sample of ethyl bromide 
treated with bromine as described above. The 
literature® reveals that ethyl bromide and ethanol 
form an azeotrope boiling at 37.6° and containing 
93 percent ethyl bromide. If our sample of ethy! 
bromide contained 7 percent ethanol, the initial, 





5M. Bodenstein and H. Liitkemeyer, Zeits. f. physik. 
Chemie 114, 208 (1924) ; M. Bodenstein and S. C. Lind, ibid. 
57, 168 (1906). 

6 International Critical Tables, Vol. III, p. 320. 
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rapid bromination would be explained quantita- 
tively if it is assumed that ethanol is brominated 
much more rapidly than ethyl bromide. In view 
of the strength of the C—OH bond (some 30 
kcal. less than that of the C—H bond), this is 
quite plausible. Probably a small amount of 
oxygen would be released and would be expected 
to inhibit bromination of ethyl bromide. There- 
fore no quantitative reliance can be placed on 
the experiments with ethyl bromide, although 
they indicate that the rate is not very different 
from that with ethane. 


CARBON-HYDROGEN BOND STRENGTH 


It is possible to utilize data presented in this 
paper to calculate with considerable accuracy the 
carbon-hydrogen bond strength in ethane. At the 
mean temperature of the experiments, 333°K, 
the heat of reaction of C2:Hs+Br—C.H;+HBr 
will be the activation energy of the reaction as 
written, diminished by the activation energy of 
the reverse, i.e., 13.3—E, kcal. While no meas- 
urements of Ey have been made, it can be esti- 
mated with close approximation. In the case of 
methane,! E, was shown to be about 1-2 kcal. 
In general it is known that such reactions occur 
on almost every collision and require very low 
energies of activation. We shall assume that in 
this case E,=0.8 kcal., whence the heat of re- 
action is 12.5 kcal. at 333°K. This has been 
recomputed to 0°K by the method described 
in I, giving AHo°=11.9 kcal. We now write: 





OF 


HYDROCARBONS 


A o° 
11.9 keal. 
85.8 kcal. 


CoH,o+ Br=C.oHs+HBr 
HBr=H+Br 


C.Hg=C2Hs+H 





97.7 kcal. 
A probable error of about +2 kcal. is assigned 


to this value. We find the following rounded 
values of this C—H bond strength 


“lk D(cotis -H) 
0 98 kcal. 
298 99 kcal. 


Steacie and Phillips’ have estimated the C—H 
bond in ethane to be 95 kcal., while Stevenson® 
has obtained the value of 96 kcal. from electron 
impact potentials. Our value of 98 kcal. would 
approach these more closely if it should develop 
that the activation energy of reaction (4) were 
larger than 1 kcal. 

Comparison with the C—H bond strength 
obtained for methane in I shows that the bond 
strength of ethane is about 3 kcal. lower, a result 
one might anticipate. It is interesting that the 
substitution of either a bromine atom or a 
methyl radical for an hydrogen atom in methane 
causes about the same decrease in carbon-hydro- 
gen bond strength. 

We wish to thank Professor G. B. Kistiakowsky 
and Dr. R. S. Halford for their interest and 
suggestions in this research. 

7E.W. R. Steacie and N. W. F. Phillips, J. Chem. Phys. 6, 


179 (1938); Can. J. Research B16, 303 (1938). 
8D. P. Stevenson, J. Chem. Phys. 10, 291 (1942). 
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An investigation has been made of the mercury (*P;) photosensitized polymerization and 
hydrogenation of butadiene. The main products of the reaction of butadiene with mercury 
(?P;) atoms are hydrogen, an acetylenic compound in the C, fraction, dimer, and polymer. In 
the presence of excess hydrogen, the products are mainly butane, butene, and octanes. The 
variation of the rate of the polymerization reaction with time, and with the pressure of 
butadiene, suggests the following mechanism: C,Hs+Hg(?P1)—~C,H.* + Heg(!So), followed by 
CyHo* + CaHo—2CuHe, or CsHo*+CaHo(CaHe)2 or CsHe*—>CyH.+ He. Polymerization arises 
largely by a free radical mechanism H2+Hg(?P1)~2H+Hg(!S9), H+CsHe C,H, CsH7 + CiHe 
—CsH:3, etc. The mechanism is shown to be consistent with the value found for the quantum 


yield. 





INTRODUCTION 


LTHOUGH the mercury [*P;] photosen- 

sitized reactions of many of the lower 
paraffins and olefins have been studied,’ no 
papers have as yet appeared on either the mer- 
cury photosensitized polymerization or hydro- 
genation of butadiene. This conjugated alkadiene 
is of very considerable interest, and the work 
described herein was undertaken in the hope 
that some light might be thrown on the mecha- 
nisms involved in its reactions. 


PART I. THE MERCURY PHOTOSENSITIZED 
POLYMERIZATION OF BUTADIENE 


Static Runs 


Both static and single-pass runs were made, 
the former at 30°C, and the latter at room tem- 
perature. 

In the earlier part of the work with a static 
system a quartz annular cell was employed 
having a volume of 342 cc. The mercury ‘“‘reso- 
nance” lamp used in conjunction with this cell 
was fashioned in the form of a trap. The trap was 
constructed entirely of Pyrex glass with the 
exception of the lower half of the outer tube, 
which was of quartz and joined to the Pyrex 
portion by means of a graded seal. The electrodes 
were sealed to the inlet and outlet tubes by 
means of uranium glass seals. The lamp was 
filled with neon at 3-mm pressure, and was 
operated at about 120 milliamp. from a 6000-volt 

*Contribution No. 1248 from The National Research 
Laboratories, Ottawa, Canada. 


1 For references, see: D. J. Le Roy and E. W. R. Steacie, 
J. Chem. Phys. 10, 676, 683 (1942). 
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sign transformer. The lamp and cell were im- 
mersed in a thermostat containing distilled 
water, and maintained at 30.0+0.1°C, the quartz 
portion of the lamp being placed in the center of 
the annular cell. The temperature was con- 
trolled in the thermostat by means of a thermo- 
regulator of the DeKhotinsky type connected to 
a knife-blade heater. It was found necessary to 
run the lamp continuously for at least 48 hours 
before its light intensity became constant. During 
this time the intensity fell to about 50 percent 
of its original value. 

In order to facilitate the determination of the 
quantum yield, and of the effect of light intensity 
on the rate, the cell system was modified to some 
extent. A cylindrical fused quartz cell, 10 cm in 
length, with windows 5 cm in diameter, was sub- 
stituted for the annular cell. Its volume was 240 
cc. The cell was immersed in a metal water bath 
provided with a fused quartz window. The dis- 
tilled water in the bath was kept at 30.00 
+0.01°C, by means of a thermoregulating system 
consisting of a mercury thermoregulator, a 
knife-blade heater, and a thermionic relay em- 
ploying a 117L7GT rectifier-beam, power ampli- 
fier tube in its circuit as recommended by 
Serfass.? The mercury resonance lamp was made 
by fusing two Pyrex-to-quartz, 13-mm O.D., 
graded seals together at the quartz ends, and 
bending the tube in the form of a U. The elec- 
trodes were fused to the Pyrex ends by uranium 
glass seals. The lamp was filled with neon at a 
pressure of 3 mm and clamped securely outside 


2 E. J. Serfass, Ind. Eng. Chem. Anal. Ed. 13, 262 (1941). 
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the thermostat about 2 cm from the fused quartz 
window. 

The butadiene used in the experiments was 
obtained from the Ohio Chemical and Manu- 
facturing Company and was shown by analysis 
to be 98.5 percent butadiene. It was further 
purified by several trap-to-trap distillations, and 
finally stored in a one-liter flask, connected to the 
manifold of the system. The storage flask was 
provided with a modified Stock valve, as designed 
by Warrick and Fugassi.* This valve proved to 
be a very convenient method of introducing 
definite quantities of gas into the cell system. 

The nitric oxide used in the work was prepared 
by the action of mercury on a 2 percent solution 
of sodium nitrite in concentrated sulphuric 
acid.! The other oxides of nitrogen present were 
removed by passing the gas through two long 
glass spirals, the first immersed in a dry-ice 
acetone mixture at —78.5°C, and the second in 
petroleum ether at — 148°C. The nitric oxide was 
then frozen out in a trap with liquid air, and 
subjected to several trap-to-trap distillations, 
the middle fraction being collected and stored 
in a similar storage flask to the one used for 
butadiene. 

Pressure changes were followed on McLeod 
gauges covering the range from 210-5 to 110 
mm of mercury. 


Analytical Methods 


The analytical methods employed underwent 
several modifications during the course of the 
work. Non-condensable gas was drawn off 
through a “pumped down” liquid air trap, by 
means of a two-stage diffusion pump, backed by 
a Toepler pump. The hydrogen in the non-con- 
densable gas was measured either by circulation 
Over copper oxide deposited on glass wool, and 
maintained at 300°C, or by transferring the 
sample to a Blacet-Leighton gas analysis ap- 
paratus, if its volume was small. 

Isolation of the fractions of various carbon 
numbers in the products was accomplished by a 
method of fractional condensation and vaporiza- 
tion, using the Ward apparatus,® as modified by 


_*E. Warrick and P. Fugassi, Ind. Eng. Chem. Anal. 
Ed. 15, 13 (1943). 
*A. Klemenc, Die Behandlung und Reindarstellung von 
Gasen (Akad. Verlagsges., m. b. H., Leipzig, 1938), p. 166. 
E. C. Ward, Ind. Eng. Chem. Anal. Ed. 10, 169 (1938). 
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Savelli, Seyfried, and Filbert.* Four traps were 
connected in series, and separated from each 
other by stopcocks. This trap system was con- 
nected at one end to the system through a 
modified Stock valve,’ and to the other end were 
connected a two-stage diffusion pump, an all- 
mercury Toepler pump, a gas measuring burette, 
and a copper oxide tube for determining hy- 
drogen. Each trap was provided with internal 
and external copper-constantan thermocouples. 
The traps were maintained at the desired tem- 
perature either by surrounding them with the 
appropriate substance at its freezing point, or 
by the modified Ward method.® In this method 
a copper cylinder wrapped with resistance wire 
was slipped over the trap, and the ends of the 
resistance wire attached to binding posts on a 
control panel. These binding posts were con- 
nected to General Radio 200-CM Variac trans- 
formers, giving a variable voltage from 0-135 
volts. The copper cylinder was then filled with 
a light gasoline. A glass cylinder, filled with 
gasoline, and a Dewar flask were then slipped 
concentrically over the copper cylinder. The 
Dewar flask was kept filled to the top with liquid 
air during the analysis. By a careful adjustment 
of the Variac it was found possible to keep the 
temperature constant to 0.5°C over the range 
from —10°C to —175°C. 

The acetylene compound found in the C, frac- 
tion was estimated by the method of Ross and 
Trumbull,’ as modified by Le Roy and Steacie* 
for their determination of acetylene. Butadiene 
did not react with the neutral 2.5 percent silver 
nitrate solution, used to absorb the acetylene 
compound. 

The temperatures used in the traps and the 
cutting pressures are given in reference 6. The 
cutting pressures were determined with a 
McLeod gauge. An atmospheric Toepler pump 
was attached to the system for transferring frac- 
tions for analysisin the Blacet-Leighton apparatus. 


Single-Pass Runs 


About a dozen single-pass runs were made in 
order to study the initial products of the reac- 


Ka," Savelli, W. D. Seyfried and B. M. Filbert, Ind. 
Eng. Chem. Anal. Ed. 13, 868 (1941). 

7W. H. Ross and H. L. Turnbull, J. Am. Chem. Soc. 
41, 1180 (1919). 

8D. J. Le Roy and E. W. R. Steacie, J. Chem. Phys. 9, 
829 (1941). 
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tion without the complications arising from 
secondary reactions. The system consisted of a 
large trap connected to a quartz annular cell 
through a fine capillary. The cell contained a 
drop of mercury, and was attached to the system 
with ground glass joints in order to facilitate 
removal for cleaning. The lamp used was a 
Hanovia Sc 2537, of the double-U type. It was 
inserted through the annulus. The cell was con- 
nected in turn to the analytical unit described 
above. 

In these runs the butadiene was first frozen 
out with liquid air in the large trap, after which, 
the system was pumped down to a pressure of 
10-5 mm of mercury, in order to remove any 
traces of air or other permanent gases that might 
have been present. The temperature of the trap 
was then raised to —78.5°C, by immersing it in 
dry-ice acetone. Butadiene has a vapor pressure 
of 11.2 mm of mercury at this temperature. The 
system was then left for six hours in order that 
the vapor pressure of the mercury might attain 
its equilibrium value. To start the run the lamp 
was turned on and the valve connecting the cell 
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Fic. 1. Pressure vs. time for a typical static run. 


to the analytical unit was opened. The butadiene 
flowed in turn through the cell, then through a 
trap at —85°C, and finally was condensed in a 
liquid air trap. The dimer formed in the reaction 
was frozen out in the trap at —85°C. The 
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acetylene compound condensed with the buta- 
diene in the liquid air trap. The mercury vapor 
pump on the analytical unit served to pump off 
the hydrogen formed, which was collected in a 
gas burette and analyzed by combustion over 
copper oxide at 300°C. 
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Fic. 2. The initial rate vs. the initial pressure for 
the static runs. 


RESULTS 
Static Runs 
The Effect of Pressure on the Rate 


Figure 1 shows a typical one-hour run at an 
initial pressure of 10 mm of butadiene. As can 
be seen from the diagram, the curve consists of 
three portions. There is a period of induction AB, 
during which time the pressure decrease is very 
small. This is followed by a linear portion BC, 
and the slope of this line is considered as the rate 
of the over-all reaction. Finally CD represents a 
falling off of the rate, caused, as has been shown, 
by the accumulation of polymer on the window of 
the cell. The length of the period of induction 
AB is directly proportional to the initial pressure 
of butadiene in the system, being about 10 
minutes at 10 mm and increasing to about 28 
minutes at 30 mm. 

Taking, then, the rate of the over-all reaction 
as the slope of the portion BC of the curve, this 
value is plotted in Fig. 2 against the initial 
pressure of the butadiene. The fact that the 
rate of the reaction as defined in this manner 
reaches a maximum at about 2 mm is un- 
doubtedly caused by the incomplete quenching 
of the Hg(*P:) atoms by the butadiene at lower 
pressures. 

That the falling off of the rate in the CD 
portion of the curve in Fig. 1 is caused entirely 
by the accumulation of a polymer on the window 
of the cell can be shown by pumping off the 
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butadiene after a run, and then making another 
run at the same pressure, without cleaning the 
cell. The rate is reduced to a value equal to the 
slope of the CD portion of the curve. 

The polymer, which was clearly discernible on 
the window of the cell through which the light 
entered, was removed between runs with a hot 
saturated solution of chromic oxide in concen- 
trated sulphuric acid. In certain cases it was 
found necessary to flame out the cell in order to 
ensure the complete removal of polymer. This 
was imperative in the case of the runs made 
with nitric oxide-butadiene mixtures, suggesting 
that the character of the polymer was different 
in this case. 

Kublitskii? found that the polymer formed 
upon illuminating butadiene sealed in quartz 
ampules with light from the mercury arc, ab- 
sorbed \2537, in agreement with our observations. 
Another interesting characteristic of the polymer 
is that it seems to exhibit the phenomenon of 
“dark growth.” If the butadiene and products 
are pumped off at the end of a run, and the cell 
again filled with butadiene, without removing 
the cell from the system for cleaning, it was found 
that the rate for the second run decreased the 
longer the cell was left before making the second 
run, 


The Effect of Nitric Oxide on the Rate 


Rate-inhibition by nitric oxide is usually taken 
to be an indication that free radical chains play 
an important part in a reaction. For this reason 
the effect of nitric oxide on the rate of the 
butadiene reaction was studied. Figure 3 shows 
a typical run with a mixture of 8 vol. percent of 
nitric oxide in butadiene. There is no period of 
induction in the reaction, and the falling off of 
the rate is caused by an accumulation of polymer 
on the wall. The rate of the reaction in Fig. 3, 
as measured by the rate of pressure decrease, is 
about 2} times as great as the rate with pure 
butadiene at the same pressure. 

The rate was measured over a concentration 
range of 0.2 percent to 20 percent of added 
nitric oxide, and it was found that the rate of 
this reaction, with a constant partial pressure 
of butadiene, was directly proportional to the 


*K. M. Kublitskii, Sintet. Kauchuk (U.S.S.R.), No. 7/8, 
30-32 (1936). 
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amount of nitric oxide present. Analysis of the 
products of the reaction showed that the nitric 
oxide had been used up in the course of the reac- 
tion. All the above facts seem to suggest that 
we are dealing here with the mercury photo- 
sensitized reaction of nitric oxide with butadiene, 
and for our purpose, therefore, these results are 
irrelevant. 


THE PRODUCTS OF THE REACTION 
(a) Static Runs 


Small amounts of hydrogen and ethane were 
formed in the reaction, together with a large 
amount of a dimer and of course the polymer 
which forms on the window of the cell, through 
which the light entered. 

Analysis of the gas non-condensable in liquid 
air from several one-hour runs, made at an initial 
pressure of 25 mm of butadiene, showed 0.06 
mole of hydrogen formed for every mole of 
butadiene decomposed. Two-hour runs at the 
same initial pressure yielded only 0.01 mole of 
hydrogen for each mole of butadiene decom- 
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Fic. 3. Pressure vs. time for a static run with 8 vol. per- 
cent of nitric oxide and 92 vol. percent of butadiene. 





posed, showing that the hydrogen formed initially 
was being used up in the course of the reaction. 

By freezing out the condensable products in 
liquid air in a small bulb attached to a capillary 
manometer, and then pumping off the butadiene 
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from dry-ice acetone at —78.5°C, a colorless 
dimer was isolated. By surrounding the bulb 
with various freezing mixtures a series of vapor 
pressure measurements were made after the 
manner described by Benson.’° Measurements 
were made from —78.5°C to 50.0°C. An extra- 
polation of the plot, log (vapor pressure) vs. 1/T 
to 760 mm gave a value of 95°C for the normal 
boiling point of the dimer. Gravimetric micro- 
analysis gave C=88.2 percent, H=11.7 percent. 
Calculated for CsH12, C=88.8 percent, H=11.2 
percent. 

The dimer showed a great tendency to poly- 
merize both in vacuum and in air. In both cases 
a yellow transluscent solid is formed, which 
slowly loses the “double bond” odor of the 
original dimer. Oxygen does not seem to affect 
the rate of the polymerization. 

About 0.6 mole of dimer were formed for every 
mole of butadiene decomposed. This was the 
average of 5 two-hour runs at a pressure of 25 
mm. 

A small amount of an acetylenic compound 
was found in the C, fraction by using the method 
of Ross and Trumbull,‘ as adapted by Le Roy 
and Steacie.* This may be vinyl acetylene. 

The amount of the polymer formed decreases 
markedly with increasing pressure. At an initial 
pressure of 100 mm there is no falling off of the 
rate even after two hours illumination. Further- 
more, there is no visible evidence of polymer on 
the window of the cell. 


(b) Single-Pass Runs 


The products from these runs were hydrogen, 
an acetylenic compound in the C, fraction, a 
dimer, and a trace of polymer. The mean of four 
12-hour runs gave: 


Products Miullimoles 
Hydrogen 2.0 
Acetylenic c.p.d. 0.5 (as C4H,) 
Dimer 0.37 (as CsH 12) 
Polymer trace. 


The Effect of Light Intensity 


The effect of light intensity on the rate was 
determined by inserting a calibrated copper 


10S. W. Benson, Ind. Eng. Chem. Anal. Ed. 14, 189 
(1942). 
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screen between the lamp and the cell, which 
reduced the light intensity to 0.32 of its full 
value. The rate at this intensity was 0.39 of the 
value obtained by using the full output of the 
lamp. From this it can be seen that the rate is 
approximately directly proportional to the light 
intensity. 


The Quantum Yield 


A rough determination of the quantum yield 
was made by comparison with the initial pressure 
rise in the Hg-photosensitized reaction of 
ethylene.’ Le Roy and Steacie found that the 
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Fic. 4. The effect of pumping off the hydrogen 
in a static run. 
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quantum yield of the ethylene reaction was 0.37 
at 13-mm pressure of ethylene.’ The quantum 
yield of the butadiene reaction calculated on this 
basis was 0.8+0.5 at 25-mm pressure. According 
to Gee, as quoted by Bolland and Melville," the 
quantum yield of the mercury photosensitized 
reaction of butadiene is in the neighborhood of 
unity. This is in agreement with our results 
within the limit of error. 


The Effect of Hydrogen on the Rate 


Two sets of runs were made in order to find 


out whether any relation existed between the 


11 J. L. Bolland and H. W. Melville, Proc. Rubber Tech. 
Conf., No. 90, p. 248 (1938). 
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period of induction and the hydrogen found in 
the products. Figure 4 shows a typical run in 
which the hydrogen has been pumped off after 
the reaction had reached the linear part of the 
curve. It is obvious from Fig. 4 that a new 
period of induction is begun when the reaction 
is resumed after pumping off the hydrogen. 

Another series of runs was made in which a 
small amount of hydrogen was added to the 
butadiene before starting the run. The period of 
induction was completely eliminated in these 
runs. The same effect is produced by stopping 
the reaction in its period of induction and 
adding a small amount of hydrogen. 

The slope of the linear portion of the curve 
increased with the percentage of added hydrogen. 
This is in accordance with the results from the 
hydrogenation of butadiene with Hg[*P; ] atoms; 
the rate of this reaction for a 10:1 ratio of 
hydrogen to butadiene is about ten times as fast 
as the reaction with pure butadiene. 


Discussion 
(a) Static Runs 


Any mechanism proposed to explain the reac- 
tions of butadiene photosensitized by Hg(*P1) 
atoms must be consistent with the following 
facts. 

(a) There is a period of induction in the reac- 
tion the length of which is directly proportional 
to the initial pressure of the butadiene. 

(b) Small amounts of hydrogen and an 
acetylenic Cy, compound are found in the 
products of the reaction. The hydrogen disap- 
pears slowly upon prolonged illumination. 

(c) If the reaction is stopped, after the linear 
portion of the pressure-time curve has been 
reached, and the hydrogen pumped off from 
liquid air, the period of induction reappears 
when the reaction is resumed. 

(d) Adding a small amount of hydrogen to 
the butadiene eliminates the period of induction, 
and the resultant rate increases with increasing 
amounts of added hydrogen. 

(e) The rate as defined above decreases slowly 
with increasing pressure. 

(f) Appreciable amounts of a liquid dimer are 
formed. 

(g) A highly unsaturated polymer is formed 
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which absorbs 42537. The amount of the polymer 
formed decreases with increasing pressure. 

(h) The quantum yield of the reaction is 
approximately unity. 

By analogy with the results found by Le Roy 
and Steacie” for the mercury photosensitized 
reactions of ethylene, the formation of an 
activated molecule scems the most probable 
primary step: 


C.Het+He(@P))>CyHe* +He (4S). (a) 


Indeed, the formation of an activated butadiene 
molecule of relatively long life by the collision of 
a butadiene molecule with a mercury (*P;) atom 
would seem even more likely than with ethylene 
since butadience has a larger number of degrees 
of freedom over which the energy imparted upon 
collision might distribute itself. 

Since the C—H bond strength in the olefins 
probably is not much different from that in the 
paraffins, the reactions 


C,He+He(?@P1)-C,H;+H+Heg('So)  (b) 
and 


C,H,*-C,H;+H (c) 


may occur to some extent, although these reac- 
tions do not seem to be of importance as initial 
steps at room temperature. 

The period of induction in the reaction is pre- 
sumably the result of competition between a 
pressure-increasing and a_ pressure-decreasing 
reaction. By analogy with the ethylene reaction, 
and in the light of our own findings, these 
would be: 

CyHe* > CyHi+He (d) 
and 
CyHe* +CyHe—(CyHe)o. (e) 


Reaction (d) would explain the small amounts 
of hydrogen found in the products, and the 
presence of an acetylenic compound in the C, 
fraction. The acetylenic compound would most 
probably be vinyl acetylene. The dimer formed 
by reaction (e) would explain how the liquid 
dimer is formed which was found in the products. 
According to Wasserman,” the Diels-Alder bi- 
molecular dimerization of butadiene results in 
the formation of 4-vinylcyclohexene. Vaughan" 
isolated what appeared to be a bicyclooctadiene. 


#2 A. Wasserman, J. Chem. Soc., p. 612 (1942). 
13W. E. Vaughan, J. Am. Chem. Soc. 54, 3863 (1932). 
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It is probable that the dimer formed in the 
photosensitized reaction is of a similar nature to 
these compounds. 

Reactions (d) and (e) are probably the only 
reactions of importance occurring during the 
period of induction which would have any effect 
upon the over-all pressure. Increasing pressures 
would favor reaction (e), and consequently the 
pressure decrease in the induction period is more 
apparent at higher pressures. The quantum 
yield of 0.8, together with the decrease in rate 
with increasing pressure, strongly suggests that 
a deactivating reaction also occurs: 


CyHo* +CyHe—2CaHe, (f) 


as in the case of the ethylene reaction (8). 
During the period of induction, then, reaction 
(d) will result in the accumulation of hydrogen 
in the system. Eventually the hydrogen will 
reach a high enough concentration to quench the 
Hg(*P;) atoms to an appreciable extent: 


H2+Hg(@P:) 2H +Hg('S»). (g) 


The hydrogen atoms formed by (g) will then be 
used up by hydrogenating the butadiene: 


CyHe +H-CyH;. (h) 


By analogy with the ethylene reaction (h) should 
have an activation energy of 5 kcal. or less. 
Polymerization would then occur by 


C,He+CsuH;—-CsHi3, etc. (i) 


The end of the induction period will occur 
when the hydrogen has built up to a high enough 
concentration to quench the Hg(§P;) atoms to 
an appreciable extent, by reaction (g). The fact 
that a very small amount of hydrogen will bring 
about a large pressure decrease, and result in 
the formation of considerably more polymer is 
additional and convincing evidence that at least 
part of the polymerization occurs by means of a 
free radical reaction. 

The above mechanism is also in accordance 
with the fact that the length of the induction 
period increases with increasing pressure. At 
higher pressure it would take a proportionately 
longer time for the hydrogen formed by (d) to 
attain a sufficiently high concentration to be 
able to quench an appreciable number of the 
Hg(*P;) atoms, and thereby initiate polymeriza- 
tion, by (g), (h), and (i). 


AND 


BE. W. KR. STSACTE 

The runs done on the effect of hydrogen upon 
the rate also afford a striking confirmation of the 
above mechanism. The fact that the period of 
induction can be reproduced by pumping off the 
hydrogen present, when the reaction has reached 
a stage where the pressure is varying linearly 
with the time, shows that the main part of the 
polymerization must be initiated by H atoms, 
as in the case with the mercury photosensitized 
ethylene reaction (8). 

Moreover, in the light of the above mechanism, 
one would expect that at higher pressures the 
polymerization would be depressed in favor of 
the dimerization. This would result in a smaller 
yield of polymer at higher pressures, a larger 
yield of dimer, and a decrease in rate. All this is 
in complete agreement with our findings. The 
quantum yield will always be small since deac- 
tivation by (f) will always remove a large pro- 
portion of activated molecules. 

The small amount of ethane identified in the 


products may have been formed by some 
sequence such as 
C,H; +H-2C2H,, (3) 
Co.H4,+H-C2H;s, (k) 
CoH3+C4H;- CoH + CaH es. (1) 


(b) Single-Pass Runs 


In the single-pass runs the rate of flow of 
butadiene through the cell was so adjusted that 
the reaction remained essentially in its period of 
induction throughout the course of the run. 
Under these conditions the main reactions would 
be (a), (d), (e), and (f). The hydrogen formed by 
(d) would never reach a high enough concen- 
tration in the cell to quench the Hg(*P;) atoms 
to any appreciable extent. Therefore one would 
expect that the ratio 


moles of dimer formed /moles of hydrogen formed 


would be small in the single-pass runs, and large 
in the static runs, where the hydrogen is being 
continuously used up by (g). Our analyses are 
in accord with this conclusion ; 


Moles dimer formed 
Moles hydrogen formed 

Static runs (1 hour) 10 
Single-pass runs 0.2. 
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Fic. 5. Pressure vs. time for a static run with 15 vol. per- 
cent of butadiene and 85 vol. percent of hydrogen 


The trace of polymer formed in the single-pass 
runs suggests that, under these conditions, poly- 
merization may have been initiated to some 
extent by C,H; radicals resulting from reactions 


(b) and (c). 


PART II. THE MERCURY-PHOTOSENSITIZED 
HYDROGENATION OF BUTADIENE 


The work previously described on the poly- 
merization of butadiene, sensitized by Hg(*P;) 
atoms, strongly suggested that the polymeriza- 
tion of butadiene was initiated by hydrogen 
atoms in the following manner: 


CsHo+He(?P1)>CaH6* + Hg('So) (1) 


CyH.*-CyHs+H, (2) 
H.+Hg(?P;) 2H + Hg ('S») (3) 
CyHe+H-CyH; (4) 
CyH;+CsHe—CsHi3— polymer. (5) 


The work described herein on the mercury 
(*P1) photosensitized hydrogenation of butadiene 
was undertaken in order to gain additional 
information on the polymerization reaction. 


Experimental 


The experiments were carried out, at 30°C, in 
a static system similar to the one employed above. 
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Fic. 6. Pressure vs. time for a static run with 9 vol. percent 
of butadiene and 91 vol. percent of hydrogen. 





The hydrogen used in the runs was prepared 
by passing hydrogen from a cylinder through 
two liquid-air traps, and then through a com- 
bustion tube containing platinized asbestos at 
500°C. It was stored in a two-liter flask provided 
with a modified Stock valve.* 

The mixtures of butadiene and hydrogen were 
made up by first introducing the butadiene into 
the cell system and measuring its pressure on a 
large McLeod gauge. The hydrogen was then 
admitted after the butadiene had been frozen 
out in liquid air. The partial pressure of hydrogen 
in the mixture was taken as the difference 
between the total pressure and the pressure of 
the butadiene. 


THE EFFECT OF TOTAL PRESSURE AND 
HYDROGEN CONCENTRATION ON 
THE OVER-ALL RATE 


Figures 5 and 6 show two runs, the first made 
with a fivefold, and the second with a tenfold, 
excess of hydrogen. There is no period of in- 
duction, and the initial rate is about ten times as 
great as the rate obtained when no hydrogen has 
been added. The falling off of the rate in the BC 
portion of the curve in Fig. 5 was shown to be 
caused by the accumulation of a polymer on the 
window of the cell through which the light 
entered. This was done by adding more butadiene 
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at the end of the run and then re-illuminating the 
mixture. The resulting rate corresponded to a 
continuation of the BC portion of the curve. 

The falling off of the rate in Fig. 6 was caused 
not by an accumulation of polymer on the 
window, but bya complete conversion of the buta- 
diene present to saturated products. Upon adding 
butadiene at the end of this run, a rate was ob- 
tained which was only slightly less than the 
initial rate. An analysis of the products of this 
run is given below. Only a few patches of polymer 
were found on the window of the cell. 

The effect of hydrogen concentration on the 
rate is shown in Fig. 7. The initial rate is here 
plotted against mm of added hydrogen, the 
partial pressure of the butadiene being main- 
tained constant at 10 mm. The relation is roughly 
linear for mixtures containing less than 80 vol. 
percent of hydrogen. At higher concentrations 
the rate falls off as would be expected. 

Increasing the total pressure at constant com- 
position did not perceptibly affect the rate up to 
a total pressure of 300 mm. However, the 
amount of polymer formed was markedly affected 
—so much so that at 300 mm the accumulation 
of polymer was sufficient to render the window 
opaque to 2537, and stop the reaction, when 
only 5 percent of the butadiene had reacted. 


The Products of the Reaction 


The products of the reaction were separated 
into fractions of the same carbon number by the 
analytical technique previously described. The 
C, fraction was analyzed by passing the gas over 
maleic anhydride at 100°C, to constant volume, 
and then determining the butene in the residue 
by bromination, using the method of Uhlig and 
Levin, in the microadaptation of Benson." 

The products of the runs made with mixtures 
containing more than 90 mole-percent of hy- 
drogen consisted mainly of butane with smaller 
amounts of methane, ethane, and octanes, 
together with a small amount of polymer. Anal- 
ysis of the run shown in Fig. 5 gave 10 percent 
methane, 1 percent ethane, a trace of propane, 
70 percent butane, 5 percent butene, 10 percent 
octanes, and a small amount of polymer. 

Qualitative analyses were carried out on the 
products of the runs from mixtures containing 
less than 90 percent hydrogen, and these showed 
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that the products were preponderately butene 
and butane, and polymer. 


Discussion 


The quenching cross section for hydrogen is 
8.6X10-'* cm?, while that for butadiene should 
not be very far from 40X10~'* cm?," by analogy 
with other unsaturated substances, although it 
has never been experimentally determined. 
Using these values in the standard kinetic theory 
expression for the collision number, we obtain 


Z(H2, Hg*)/Z(CsHs, Hg*) = [H2]/LCsHe J, 


where [H2] and [C,He¢] are the molar concen- 
trations of hydrogen and butadiene, respectively. 
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Fic. 7. The intial rate vs. mm of hydrogen added to 
10 mm of butadiene. 


In other words, the extent of quenching of 
Hg(*P1) atoms done by butadiene and hydrogen 
in a mixture will be approximately proportional 
to the ratio of the molar concentration of the two 
substances in the mixture. It follows, then, that 
in a mixture of butadiene and hydrogen con- 
taining 90 mole-percent of hydrogen, the hy- 
drogen will do about 90 percent of the quenching. 

Bearing the above facts in mind, two primary 
steps seem likely; namely, 


H.+Heg(*P:)—2H +Heg(!Sp), (1 ) 


C,H. +Heg(?P1)-C.He*. (2) 
4A. C. G. Mitchell and M. W. Zemansky, Resonance 


Radiation and Excited Atoms (Cambridge University 
Press, 1934), p. 224. 
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The large amounts of polymer formed in the 
hydrogenation runs lend weight to our previ- 
ous conclusion that the polymerization occurs 
through the addition of an H atom— 


CysHe +H—-C,H,;, (3) 
CyHe+CsH;-CsH1;—polymer. (4) 


Another step which is energetically possible and 
may be important is, 


C,H6*+H2-C,H; +H. (5) 


Certainly the fact that the main product formed 
is butane in a mixture containing 90 mole-percent 
of hydrogen would seem to indicate that there is 
a relatively high concentration of both C,H; 
radicals and H atoms. The butane formed, under 
these conditions, would likely arise from some 
such sequence as 


C,H; +H-C,Hsg, (6) 
C,Hs+H—-CyHyg, (7) 
followed by 
C4Ho+CyH;-C,He+CuH w, (8) 
and/or 
C,Ho+He-CyHwt+H. (9) 


Disproportionation is not considered impor- 
tant with normal radicals; however, the type of 
cross disproportionation represented in (8) may 
be quite important, owing perhaps to the C—H 
bond in the C,H; being slightly weaker than it 
is in C4Ho. 

Reactions of radicals with hydrogen molecules 
seem to possess rather high activation energies, 
probably of the order of 8 to 15 kcal., and 
therefore one would expect (9) to be relatively 
slow. However, a butyl radical formed by the 
successive addition of three H atoms to butadiene 
might conceivably have been able to retain some 
of this energy of hydrogenation, and, for this 
reason, we might be dealing in (9) with an 
activated radical. Of course, if butene is formed 
as an intermediate, and this seems to be the case, 
it would quench almost as strongly as butadiene, 
and consequently direct photosensitized hydro- 
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genation might occur. The octanes formed un- 
doubtedly arise from radical recombination of 
butyl radicals. It is of interest to note here that 
Moore and Taylor’ found that octanes were the 
main products in the mercury photosensitized 
hydrogenation of n-butene, while we obtained 
only a 10 percent conversion of butadiene to 
octanes, our yield of butane being correspond- 
ingly greater. This discrepancy could be ex- 
plained if we assume that a large percentage of 
the butyl radicals formed from butadiene by 
hydrogenation remained in an activated state 
long enough to collide with a hydrogen molecule, 
forming butane as in (9). 

Presumably the methane and ethane formed 
could be accounted for by a series of “‘atomic 
cracking’”’ reactions starting with the butyl 
radical. However, a detailed study of these 
secondary reactions is beyond the scope of the 
present investigation. 


Conclusions 


The present investigation on the mercury 
photosensitized reactions of butadiene seems to 
suggest the following conclusions: 

(a) Collisions of butadiene molecules with 
Hg(*P;) atoms result in the formation of ac- 
tivated molecules of relatively long life. 

(b) These activated molecules are mainly used 
up by the following reactions: 


C,He*-CyHit+Hs, (i) 
CyHo*+CaHo—2CaHe, (ii) 
CyHoe* + CaHo—(CaHo)2. (iii) 


(c) Hydrogen atoms are formed when the hy- 
drogen formed by (i) reached a high enough con- 
centration in the system to quench the Hg(*P;) 
atoms appreciably. 

(d) Polymerization then occurs by the reac- 
tions: 


CyHe+H—-CyH;, 
CyHet+ C,H;-CsHi3, etc. 


(iv) 
(v) 








1 W. J. Moore and H. S. Taylor, J. Chem. Phys. 8, 504 
(1940). 
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Appropriate coordinates are set up from the standpoint of group theory for describing the 
normal modes of oscillation in such a manner that maximum factorization of the secular deter- 
minant is accomplished. The cubic and quartic portions of the anharmonic potential function 
are derived. The components of vibrational angular momentum are set down. The complete 
valence-type potential function is discussed; explicit relations are derived between the gen- 
eralized force constants occurring in the secular determinant and the valence force constants 


for CHeDs. 


1. INTRODUCTION 


HE tetrahedral X2YZ» molecular model has 
symmetry C2, and includes, as examples, 
the dihalogen substituted methanes, CHF», 
etc., the dideutero substituted methane, CH2Dsz, 
SO.Cle, and others. The normal vibrations of 

















this molecular model and the quantum-mechan- 
ical evaluation of its allowed rotation-vibration 
energies have not been treated extensively in the 
literature. Rosenthal'? has studied the calcula- 
tion of the normal frequencies of CHF, and 


~1J. E. Rosenthal, Phys. Rev. 45, 538 (1934). 
* J. E. Rosenthal, Phys. Rev. 46, 730 (1934). 
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CHD». type molecules, Dennison* has discussed 
the frequencies of CH2De, and Ta-You Wu? has 
given a qualitative description of the normal 
modes. References to experimental papers on the 
Raman and infra-red spectra of tetrahedral 
X:YZ2 molecules are given by Dennison and 
Ta-You Wu.*4 

This paper is intended to form the basis for 
the interpretation of both the gross vibrational 
and fine rotational structure of the Raman and 
infra-red spectra of tetrahedral X2YZ»2 mole- 
cules. Part I deals with the purely vibrational 
aspects of the problem and includes discussions 
of the molecular symmetry, characteristic modes 
of vibrations and potential functions. Part I] 
deals with the evaluation to second-order ap- 
proximation of the allowed energies of the rotat- 
ing vibrating model and includes discussions of 
anharmonicities, coriolis interactions between 
rotation and vibration, dependence of the rota- 
tional constants on the state of vibration, effects 
of centrifugal stretching, etc. 


2. GEOMETRY AND SYMMETRY 


The equilibrium configuration of the model 
assumed for the tetrahedral XYZ» molecule is 
shown in Fig. 1. A right-handed body-fixed 
rectangular coordinate system xyz, with origin 

3D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 


*Ta-You Wu, Vibrational Spectra and Structure of 
Molecules (Prentice-Hall, Inc., New York, 1939), p. 244 ff. 
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TABLE I. 
atom r mr Xor Yor Zor 
X 1 m ao 0 bo 
X 2 m — do 0 bo 
Z 3 n 0 aon do —€o 
Z 4 n 0 do —<€o 
y 5 M 0 0 Co 








at the center of mass of the molecule, is adopted 
for specifying the equilibrium positions of the 5 
atoms. The coordinate system is so oriented that 
its axes are the principal axes of inertia. The 
equilibrium coordinates and masses of the re- 
spective atoms are given in Table I. 

The model belongs, in the equilibrium con- 
figuration, to the symmetry point group C2, 
whose characters are given in Table Il. The 
symmetry elements of the point group shown in 
Fig. 1 are a twofold axis of rotation and two 
mutually perpendicular reflection planes inter- 
secting in that axis. The symbol E denotes the 
identity operation, C, rotation by m7 around the 
twofold axis, and o and o’ denote reflection in 
the respective planes. The symmetry types A 
and B denote quantities which are, respectively, 
symmetrical and antisymmetrical under rotation 
about the axis C2; subscripts 1 and 2 denote 
quantities which are, respectively, symmetrical 
and antisymmetrical under reflection in the 
o-plane. Several authors®~® have described meth- 
ods of obtaining many properties of a molecular 
model from the table of symmetry characters. 
All the states of the XYZ,» model are non-degen- 
erate. The total angular momentum component 
P, is of type Bs, P, of By, and P, of A». The 
coordinate x belongs to type Bi, y to Bs, and z to 
A,. There are nine normal modes, 4 of which 
belong to type A;, 1 to A», 2 to Bi, and 2 to Bo. 
The components of the dipole moment behave 
like x, y, and z; hence, only vibrations of types 
A,, By, and Be will have fundamentals active in 
infra-red. Any combination or overtone of type 
A, will be inactive in the infra-red. All classes 
of vibrations have fundamentals, overtones, and 
_>E. Wigner, Gottinger Nachrichten, p. 133 (1930); 
Gruppentheorie (Braunschweig, 1931). 

»L. Tisza, Zeits. f. Physik 82, 48 (1933). 

: E. B. Wilson, Jr., J. Chem. Phys. 2, 432 (1934). 

J. E. Rosenthal and G. M. Murphy, Rev. Mod. Phys. 
8, 317 (1936). 


* A. G, Meister, F. F. Cleveland, and M. J. Murray, Am. 
J. Phys. 11, 239 (1943). 


combinations which are Raman active; vibra- 
tions belonging to A; are polarized and the 
others are depolarized. 


3. INTERMEDIATE VIBRATIONAL COORDINATES 


The instantaneous positions of the atoms are 
described during vibration by a right-handed 
body-fixed rectangular coordinate system, which 
coincides in the equilibrium configuration with 
the system described above. The instantaneous 
position of the rth atom is given by 


ce. Vay Zr) — (xor+x,, Yor+yrs Zor+2,), 


where x,’, y,’, and z,’ denote components of dis- 
placement from the equilibrium position. There 
are 3N(N=number of atoms) displacement 
coordinates x,’, y,’, and z,’ of which only 3N—6 
=9, corresponding to the number of degrees of 
vibrational freedom, are linearly independent. 
Six linear relations among the displacement 
coordinates are given by the Eckart condi- 
tions:'°" (a) the center of mass remains at the 
origin during vibration, >> m,x,’=0, etc., and 
(6b) the angular momentum of vibration vanishes 
in zero-order approximation, 


Poz = Dm, (YorZr’ —ZorVr') =0, ete. 
For the tetrahedral X:YZ. model the Eckart 
conditions are: 
m(x1' +x") +n(x3' +x4') + Mx,’ =0, 
m(y1' +2") +n(ys' +94’) +Mys'=0, (la) 
m(2y' +22") +n(23' +24’) + Mz,’ =0, 











TABLE II, 

E C2 o o’ 
ioe 1 1 1 
42] 1 1 -1 —1 
By} 1 —1 1 —1 
Bz} 1 —1 -1 1 











Por = ndy(s4' — 2’) — mbo(y1' +2") 
+neo(ys' +s) — Mcoys’ =0, 
Poy = mbo(x1' +.X2') — neo(x3' +24) (1b) 
+ Mcoxs’ +mao(2' — 2’) =0, 
02 = Mao(y1' — yo’) + ndo(x3’ —x4/) =0. 
1 C, Eckart, Phys. Rev. 47, 552 (1932). 


1 E. B. Wilson and J. B. Howard, J. Chem. Phys. 4, 
260 (1936). 
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A suitable set of 9 linearly independent inter- 
mediate generalized coordinates, which are 
mutually orthogonal and of proper symmetry, is: 


£1 = (21'+22')/2, 

£2 = (23' +24’) /2, 

£3 = (x1' — x2’) /2, 

£,= (ys' — ya’) /2, 

05 = (x3 — x4’) /2, (2) 

no = (x1' +x2')/2, 

m7 = (x3' +24’) /2, 

63= (n' +32")/2, 

o= (ys +4’) /2, 
where £ denotes a coordinate belonging to sym- 
metry type A1, o to Ao, n to Bi, and ¢ to By. It 
follows from Eqs. (1a), (1b), and (2) that the 15 


displacement coordinates are the following func- 
tions of the generalized coordinates: 


x1’ = net és, 

Xe = 6 — £3, 

X3’=m+45, 

x4’ = 17 — 95, 

x5’ = —(2/M)[mnet+nnz |, 

yi’ = {s— d505, 

yo = Sst os05, 

ya’ = Sot bs, (3) 

ya’ =fo— bs, 

ys’ = —(2/M)[mis+nio], 

ar’ = £1 +- dens — orm, 

20’ = £1— dens +orn7, 

we! = be— Gate t dnt 

34’ = bo+dsfs— dots, 

5’ = —(2/M)[mé,+né], 
$5 = (ndo/mao), 
$6= (bo—Co)/do, 
7 = (n/m) (e9+ 0) /ao, 
os = (m/n)(bo— Co) /do, 
bo = (€0+¢0) /do. 

The kinetic energy of vibration, 


T,=(1/2) ©, mt? +y2+ 2,7) 


ta 


where 


(3a) 
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must be of symmetry type A; and is given in 
terms of the generalized coordinates by: 


27. =m é,° + rks? + usé3” + sé + 2ur2b io 
+ sos? + mene’ + wr 7 + 2u6r Hen 
+ usts?+ pots + 2usof sto, (4) 


where 
pi =(2m/M)(M+2m), 


2 = (2n/M)(M+2n), 
u3= 2m, 
bie = (4mn/M), 
bs = (2n/ may?) (mag?+nd,?), 
Me = wi + 2mge’, 
Mz = wo +2m¢,’, 
Her = M12 — 2Mgod7, 
Ms = Mi t+2nd¢s’, 
Mo = Mo + 2ngo", 


so = Miz — 2Ndsdy. 


bg =2n, 


(4a) 


The harmonic (quadratic) generalized poten- 
tial function Up is of symmetry type A; and is 
given by: 


2 Uo= hits? + oko? + aks? + aks? + 2Riokike 
+ 2Riséiés+ 2Riséigst+ 2hostots+ 2koskoks 
+ 2Rsstsks t+ Rsos’+ Rone’ + krn7” 
+ 2Rernenz t+ kets? + Role? + 2Reofsis, (5) 


in which the k’s are generalized force constants 
which are given in Sections 8 and 9 as functions 
of the special force constants in the valence-force 
potential function. 


4. NORMAL COORDINATES 
The Lagrange secular determinant, 
\AT,— Uy! =0, 


for evaluation of the normal or characteristic 
frequencies of the model is of ninth order but 
factors into four steps, each of which is associated 
with one symmetry type. The roots \ of the 
determinant are related to the normal frequen- 
cies w (measured in cm) by the expression 
\=4r°c*w*, if c is the velocity of light in cm/sec. 
The subscripts 7, 7, &, and / are used in this paper 
in the following manner to denote the normal 
frequencies belonging to each symmetry type: 
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i=1, 2, 3, 4, for type Ai; j=5 for type Ao; 
k=6, 7 for type Bi; and /=8, 9 for type Bz. The 
subscript s is used to denote the intermediate 
coordinates, defined in Eq. (2), as follows: 
s=1, 2, 3, 4, for the é’s of type A1; s=5 for o of 
As; s=6, 7 for the n’s of By; and s=8, 9, for 
the ¢’s of Bo. 

The steps of the secular determinant associated 
with symmetry types Ai, Ae, Bi, and Bz, re- 
spectively, are as follows: 



























































Aui— ki Api — Riz — ks —ki 
Auwi2—Ri2| Auo—ke | —kos | —Res 
=0, (6a) 
—kys —kes |Aus—ks| —Rss 
—kiy —ko —ks Aus —k, 
Abs eats ks = 0, (6b) 
Aus — Re Auer — Rez 
=0, (6c) 
Auer — Ror | Aur —k 
Aus — ks Aus9 ac Reg 
=(. (6d) 
Auso— Reg | Auo— Ro 














In principle, the generalized force constants 
occurring in Eq. (4) can be evaluated if one sub- 
stitutes into the secular determinant the values 
of \ corresponding to the experimentally ob- 
served normal frequencies w corrected for 
anharmonicity, etc. In this particular case there 
are 17 generalized force constants occurring in 
Eq. (5) to be determined from only 9 normal 
frequencies. It is, therefore, necessary either that 
additional data be obtained from isotopic mole- 
cules or that the 17 generalized force constants 
be expressed in terms of not more than 9 inde- 
pendent force constants in a physical potential 
function such as one of the valence or central 
force type. The number of independent physical 
force constants can be reduced by (1) arbitrary 





omission from the physical potential function 
of certain terms assumed to be small, or (2) 
assumption that the values of certain physical 
force constants, such as the CH force constant, 
are known from previous investigations. 

The transformations from the intermediate 
generalized coordinates defined in Eq. (2) to the 
normal coordinates Q are given by: 


£.=D;.0;...s=1, 2, 3, 4;1=1, 2, 3,4, (7) 
O,=,jQ0;...... s=j=5, (8) 
Ne= >t MQ. . .S=6,7;k=6, 7, (9) 
f=): 1Q:.. .s=8, 9;1=8, 9. (10) 


The transformation coefficient 1,; is given by: 
Ns 57 = N55 = bs. (8a) 
The coefficients 1, are given by: 


Nox = (Auoz — kz) / Nz, 


(9a) 
Nz = — (Auer — Rez) / Nx, 
where 
Nu=[(Aumer — Rez)? + (Anu — Rr)? ], 
if k=6, 7. The coefficients n,; are given by: 
Msi = (Aiwa — ko) / Ni, 
(10a) 
91 = (Arps9 — Reo) / Ni, 
where 
Ni=[(Arpso — Reo)? + (Arms — Ro)? }', 
if 1=8, 9. The coefficients n,; are given by: 
Nix=Asi/ Ni, (7a) 


where A,; is the cofactor of the sth element in 
the first column of the determinant given in Eq. 
(6a) with \ set equal to Aj, 


N= Cu1A 12+ u2A 02+ ys 3:?+ mA s?+ 2u2A 1:Aoi }? 


and s=1, 2, 3, 4; <=1, 2, 3, 4. 
In terms of- the normal coordinates Q, the 
quantities T, and U) are: 


9 
T,=(1/2) © Q,, (11a) 
h=1 
and . 


9 
Up= (1/2) & AnQi?=2 7c? DS wi 7Q,?. 


h=1 h=1 


(11b) 


The approximate forms of the normal modes are 
illustrated in Fig. 2. Similar pictures of the 
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modes are given by Wu.‘ The subscripts on the 
w’s in Fig. 2 follow the notation described above 
for modes associated with the various symmetry 
types. 

The subsequent discussion can be simplified 
by the introduction of dimensionless normal 
coordinates gq, defined by the relation: 


On=(h/2mcwn)*gn, (12) 


where fA is Planck’s constant divided by 2x, and 
c is the velocity of light. In terms of q, the zero- 
order vibrational energy is given in the Hamil- 
tonian form by: 





H=(he/2) & wil (pi2/h®) +7], (13) 


h=1 





Wa (Ay) t 
where p,=0T/dg,=momentum conjugate to gy. | 


5. ANHARMONIC POTENTIAL FUNCTION 


Each term in the expression for the total 
energy must belong to symmetry type A. It can 
be shown by direct multiplication of the charac- 
ters in Table II that the cubic combinations of | 
normal coordinates occurring in the anharmonic 
portion of the potential function of a molecule 
belonging to symmetry group C., can be A;’, 
AA, B;°A,, BZA te and AoB,Bo. Likewise the . 
quartic combinations occurring are A;‘, Ao}, 
By‘, B2', AvAs’, APB, AB’, A’B,’, AB’, 
BY/B-?, and AjA 2B, Bo. 

The most general cubic portion U; of the an- 
harmonic potential function of the tetrahedral 
X:YZ»2 model is given by: 


U, he= Biget L > Biv G?qi' 


+X Xu 2 Biv in QiQirir’ 


t 











+D DL Bigg? +X LD Biegige 
: z é & 


+22 oe Bink GiQedee +L Zz Bingigr 
&* l l 


‘ I 


+2 . e Burg igiqu 
i l 
+2 DD Biergiqngs. (14) 
i ® 2 


The quartic portion U2 of the anharmonic poten- 
tial function of the tetrahedral X:YZ» model, 
which contributes to the allowed energies to Fic. 2. 
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second-order approximation,” is given by: 


Us he= ) a ViiiiGi 445° > Vid? Qi"* 
HD VissidHt+D Vivwege' 
j k 
HY DY yranrergege+Y yuugi 
& & l 


$¥0 DY yuvva’?a’ +d DY viisq2? 
l i 1 


j 


+>. > Vike? ge+ > pk Viing?gr 
i k i l 
+E DY vimeg?a?+D DL vingirae 
ik iol 
+0 Dd yecngi?g?. (15) 
. 


In Eqs. (14) and (15) the subscripts 7, j, k, / 
denote coordinates belonging, respectively, to 
symmetry types Ai, As, Bi, and By as specified 
in Section 4. Primed and unprimed subscripts in 
a given summation are used to denote different 
coordinates of the same symmetry type; for 
example, Bj; may denote 8123, 8124, Biss, and 
8234. Each combination of subscripts may occur 
only once, i.e., permutations of a given com- 
bination of subscripts are not allowed; for ex- 
ample, B123 is the same as {132, B213, etc. The coef- 
ficients 8 and y occurring in Eqs. (13) and (14) 
are expressed in cm! 


6. VIBRATIONAL ANGULAR MOMENTUM 


An internal angular momentum of vibration 
may arise when two or more normal modes of 
vibration occur simultaneously. This angular 
momentum is ordinarily of small magnitude 
unless the two frequencies involved are nearly or 
exactly equal. This type of angular momentum 
is of importance in connection with coriolis inter- 
actions between rotations and vibrations dis- 
cussed in Part I]. The components p,, p,, and 
bp: of internal angular momentum are of the 
same symmetry types as the components of total 
angular momentum and belong, respectively, to 
types Bo, » Br and A». These components are 





* The most general anharmonic potential function con- 
tains the following additional quartic terms which are of 
odd degree in one or more coordinates and do not, there- 
lore, ordinarily contribute to the allowed energies: qi*gi’, 
qi gir qi, i Gir’, QRQk'; qu°qu » Qi - qi°, Gi ae gx, Gigi "Gigk’, 
qi?Qiger, 4 iy QiQirGinr, G7QIGr, GP QK e's PQQ, Gequrgr’, 
Teqeegiqur, Qeqiqur, and gig udu. 
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defined as: . 
= 2 Yr ) ’ 


7. 7 
—X; 2; ) 


pze= Ler m,(yy' 2,’ 
py= Ler m,-(B,' a,’ 

= Le m(x,' Yr’ = He Be) 
where, ~; may contain combinations of types 
AB, and A,B2; p, may contain A2B, and A,A,, 
and p, may contain A,;A,2 and B,B,. The com- 
ponents are expressed in terms of the dimension- 


less normal coordinates gq, and _ conjugate 
momenta p, as follows: 


be = Dir Exsl(ws/wr)*geps — (we/ws)*QsPx ] 
+31 Di Eni (wr/wi)*gipi— (wi/wi)*qipi), 

by = Xi nisl (ws/1)'qips — (wi/ws)* G57] 
FXe Di mil (we/wi)*gipe— (wi/on)'gupi), 

b2= Li FL (ws/wi)*Gids— (wi/ws)*gspi ] 
+e Dei Sil (wr, wr)*gupe — (wi/ wr) *Qupr). 


The coriolis parameters or moduli of internal 
angular momentum are defined as follows: 


£15 = (2ndo/ mao?) N55 m (bo — Co) Nex — 2(€o+ Co) N71. |, 
£1; = { (urttsr t+ pier) M1 + (moor + wieMsi) M2: 

+ (2/do)[m (bo — co) ns: — 2(€o-+€0) M91 |714:} , 
nis = (2//do)[n(€o+ Co) 91 — m (bo — Co) N81 ]M55, 
mi= | —(uimert mreMz) 1; — (Meme + mieMex) No: 

+ (2/ao)[m (bo — Co) M6. — n(€o+Co)N7% |nN3:}, 
Cis= —[(2ndo/ao)n3it+ mata: Mss, 
Cer = | (uiment orem) M81 + (mets, + wien) N91} , 


(16) 


(17) 


where 7 assumes the values 1, 2, 3, 4; k the values 
6, 7; and / the values 8, 9. 


7. GENERAL VALENCE POTENTIAL 
FUNCTION 


A harmonic potential function suitable for the 
interpretation of experimental data and the cor- 
relation of vibration spectra of related molecules 
should contain force constants possessing phys- 
ical significance rather than the generalized alge- 
braic constants occurring in Eq. (5). An appro- 
priate physical potential function is the valence 
type in which the coordinates used are the 
valence bonds and their included angles. The 
valence coordinates suitable for description of 
infinitesimal vibrations of a group of three atoms 
(m,n, and a central atom C) are illustrated in 
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Fig. 3, where fm=fomt+6fm, Tn=TontOrn, and 
Bmn=Bomnt+6Bmn; Tm; Yn, ANd Bmn denote, respec- 
tively, the instantaneous values of the two 
valence bond lengths and the included bond 
angle, while rom, fon, aNd Bomn denote their equi- 
librium values and 6rm, 6%n, and 658n, denote 
changes occurring during vibration. The quantity 
58mn is given for infinitesimal vibrations by 


5Bmn _ [PomndPmn a (Tom—Ton cos Bomn) orm (18) 
= (Ton —Tom COS Bomn) 6? n_1/(YomYon sin Bomn) 


where Pmn is the distance between atoms m and n. 
In the tetrahedral X.YZ2 molecule the relative 
displacements of the various pairs of atoms 
during infinitesimal vibrations can be obtained 
in terms of the intermediate coordinates by the 
substitution of Eqs. (3) into the expression 


5P mn = L(Xom—Xon) (Xm' — Xn’) 
+ (Yom — Yon) (¥m’ — Yn’) 
+ (om — Zon) (2m! —2n') ]/Pomn. 
It follows that: 
dpi2=283, Spa= —2ks, 
dpis=Lhi—fetSstSfal/ pois, 
dpu=Lfhi—fe—fs—f J/ Pos, (19) 
bpu=Lhithetfs—fa]/ pois, 
5pes=(fitfo—fstSs]/ pois, 


where 

fi=(b0+ 0) (£1 — &) +a0&3 —doks, 

f2= (do+dods)as, 

fs=[ao+ (bo+0)¢6 ne — [aot (bo+€0)¢7 Jn7, 

fa=[do+ (bo+€0) G3 ]&s — [dot (bo+€0)o0 ]is, 

Pos = Lao? +do?+ (bo+ eo)? }?, 
in which 95, etc. are the quantities defined in 
Eq. (4a); 
§pis=6r1=([gitg2]/ro1, 


(20) 
5pes = bre = [g1— ge |/ro1, 


where 

81 = (bo— Co) (ure t+ micée) /2m ]+aoés, 

2 = dol (ui/2m) +6 ]ne+aol(u12/2m) — beb7 ]n7, 
Yor = Yoo = Pos = Lao? + (bo — co”) }*; 


and 
5pss = 6r3= [gs — £4 ]/Tos, 


21 
Spas = 5r4= [gst g4 ]/r0s, _ 
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where 
&3= — (€otco)L (mice + megs) /2n.]—doks, 
84= dol (u12/2n) — soo ]istdol (us 20) +¢s? ]i&o, 
ros = Toa = Poss = Ldo?+ (€0+ 0)? }}. 


Explicit expressions for the 68m, in terms of the 
intermediate coordinates can be obtained by 
substitution of Eqs. (19), (20), and (21) into 
Eq. (18). 

The most general harmonic potential function, 
which can be set down in terms of the valence 
coordinates, contains all possible quadratic 
combinations of these coordinates. In the so- 
called ‘‘valence’’ potential function it is ordi- 
narily assumed a good approximation to include 
squared terms in the 6r’s and 68’s and to omit 
all or part of the cross products. It has been 
found in many cases that such an approximation 
gives a fairly good interpretation of the experi- 
mental data; however, the most complete poten- 


c a 


Tn — Baan 


mn 








Cc 
Fic. 3. 


tial function possible should be used to justify 
analyses of the excellent vibrational data avail- 
able in many cases. In simple cases, as for ex- 
ample that of H,O,'*!* where evaluation of the 
constants in a complete potential function has 
been possible, it has been found that none of 
these constants is negligible. A discussion of the 
complete valence-type potential function for the 
tetrahedral XYZ: model is given below in such 
a form that approximations can be made simply 
by dropping the appropriate terms. 

The valence coordinates 68m, and 67m are not 
suitable for the description of normal modes 
because (1) they do not generally fall directly 
into the required symmetry types and (2) they 
are usually redundant in the sense that they are 
more numerous than the degrees of vibrational 

13B. T. Darling and D. M. Dennison, Phys. Rev. 57, 
128 (1940). 


14.W. H. Shaffer and R. R. Newton, J. Chem. Phys. 10, 
405 (1942). 
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freedom. It is possible, however, to set up linear 
combinations of the valence coordinates which 
do possess the appropriate symmetry properties 
and a physical potential function based on such 
coordinates can be used in the interpretation of 
data if relations are found between the valence 
force constants occurring in such a function and 
the generalized constants occurring in Eq. (5). 
For the tetrahedral X2YZ» model there are ten 
possible valence coordinates, namely, four 6r’s 
and six 68’s, while there are only nine degrees of 
freedom. The work of obtaining the expressions 
relating the k’s of Eq. (5) to the valence force 
constants, which is a formidable task if done by 
direct substitution, is simplified by the adoption 
of the following ten combinations of valence 
coordinates : 


Ay = (671; +672) /2, 

Ao = (6r3+6r4)/2, 

As = (101/2)6812, 

As= (13/2) 6831, 

Ao = (ro1%03) (6813 +5821 +6814 + 6803) /2, 
As = (101% 03)'(6813-+ 6801 — B14 — 5823) /4, 
Ae = (67; — Gr) /2, 

Az = (701703) (6813 — 5821 +6814 — 6Bo3) 4, 
As = (673— 674) /2, 

Ag = (701703)! (5813 — 5824 — 6814+ 5Be3) /4. 


The variables A;, Ao, A3, and A, are of sym- 
metry type A1; Ao, of type Aj, is taken to be the 
redundant combination and can be expressed 
uniquely as a linear combination of the other 
four variables of type Ay. As is of type A3; 
As and A; are of type B;; and As and Ag are of 
type B;. The above symmetry combinations can 
be expressed in terms of the intermediate coor- 
dinates as follows: 


Ai =A.i:4+-Biiot+Ciis+Dit; (i=1, 2,3,4,0) (23) 
where 

Ay cos an2= — A3/sin any =p, 2m, 
— As/ Cos ao33= Ay/Sin aogs = wy2/2n, 


B, ‘cos ane = — B; ‘sin Gore = fye/ 2m, 


— Bs/cos ao35= By/sin ao3y= pe/2n, 


C.=C,=0, 


Ci =sin ao12, C3= COs ao, 
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D,=D3;=0, De= —sin az, 
Ay =[ (bo +0) +08 Bors(7o013A 1+ 701A 2) 

— (79:Ai+73A2) }/do. 
Bo=[ — (bo +60) +08 Bors(70381 — 71 Be) 

— (7B, —ro3B3) | 60, 


Co= [sin 912(%93 COS Boi3 — Yor) +a | do, 


dD; = — COS o34, 
(23a) 


Do = [— sin ato3a (Yor COS Bors mae ro3) —dy | Ov, 


if Aomn = | Ee » 8 60 = (Yo1%o3)? sin B13; 


A; =Aso5, (24) 

where 
As= _ (us 2n) (ao, 60) ; (24a) 
Me=Aonet+Bon;, A:;=Asne+Bin;, (25) 

where 


(Ag /me) = (Be/ uz) = (Sin aye / 207), 
A; = [a+ 6(b0+ eo) 
+ (ps /2m)(ro3 COS Borz— 71) SIN ame | 60. 
Bz =[ — ao— $7(bo+ 0) 
+ (1467/2072) (ro3 COS Bors — 71) SIN aor2 |, 60: 
As=AststBafo, Av=Aoist Bolo, (26) 
where 
(Ag/pso) = (Bs/ us) = —sin aos; 27, 
Ag=[do+os(bo+eo) + (ps9, 277) 
X (703 — Yor COS Borz) SIN ac34 | 60, 
By =[ —do— d9(bo+ €0) + (us, 27) 


X (vo3—- ¥o, COS Bois) sin ao3s | Oo. 


(26a) 


The most general harmonic potential function 
U>y is given in terms of the symmetry com- 
binations by: 

U= Ky, Ar+ KoAe?+ K3A;°+ AyAe+ Ko-Ao" 

+ 2K yAj;Ao+ 2A 134143 4+ 2A yA Ay 

+ 2Ko3A,A3+ 2K 2;A2A\4+ 2A gy AA, 

+ 2Ko1A9Ai+ 2A 024042 + 2 AozAoA3 

+2KoAoAyt+ KsAs?+ Ag Ac? + A; A," 

+2Ke7AsA;+ KgAs? 

+ KoA? + 2K gg Aso, (27) 

in which cross products occur only among vari- 
ables of the same symmetry type. The complete 


valence-type potential function U» for the tetra- 
hedral X.YZ>» model can be written in terms of 
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2Uo= cal (6r1) + (62) | +er'[(6rs) + (671) | +-corox(6812) + ¢2'ros(6831) +2" roxreal (581s) + (58.) 
+ (68; Ry aa (5823) | + 63671672 + €3'6r3674 + 3" (671 + 6r2) (673 +574) + C4701 (671 +672) 5812 
+-c4'rea(5rs +573) 5851+ rovres {€4’’ [671 (6813+ 6814) + 5r2(5823 +6821) | +04’ [r3(6813 + 6828) 
+ 674(6811 +6604) |} + c5703(671 + 672)5835 + 65’ 701(673 +574) 5812 
+ reves és’ [dri (5824+ 5823) + 672(681s+ 6814) |t+¢5'”"[6r3(5814 + 5821) + 674(5813+ 5828) J} 
+ Poros | 653125834 + C6’ (58136821 +48116823) + C7 (58136814 + 58236804) + C7” (68136823 + 68145821) } 
+ ruvresl Cr!" ror8812-+6;'”"rox88e4 (8819+ 682s +6814 +6623) (28) 


where 

Ko=2(202"’ +66 te; te’); Ki=2eates: 
Ke =2e)' +c,’ ; K3=40; Ay=40'; 

Kg =2(202"+¢¢6 —G;—¢;'); Ke=2e1.—¢3; 


Kz = 2(2¢2" — ce +6;—«7'); 
Ko = 2(2¢2" — ce’ — +07’) ; 


Kg= 2c,’ —¢s3' ; 


Ka= 2(c4'’ +05"") : 


- wr tN) 4 ” > mr (29) 
Koo = 204’ +05"); Kog=4e7""; Kog=4e;'"; 
| Cre ye , er — 2be: fisa = Zee* 
Kos = 2c’; Ku= 2c,’ : A34= 266; 


Kez = 2(c4/' — cs’) ; Agg= 2(c4/"" = “>. 


Substitution of Eqs. (23)—(26) into Eq. (27) and 
comparison of the resulting expression with Eq. 
(5) lead to the relations between the K’s of 
Eq. (27) and the generalized force constants 
k, of Eq. (5); substitution of Eq. (29) into the 
latter relations yields the expressions for the 
k.,, in terms of the valence force constants. 


8. VALENCE POTENTIAL FUNCTION FOR 
CH;D., 

The CH2D» molecule is the special case of the 
tetrahedral XYZ. model for which, in the 
equilibrium configuration, all the valence bonds 
have the same length 7p as in CH,y, and all the 
bond angles have the tetrahedral value Bo, 
=109°27’. It follows for CH2D, that the equi- 
librium dimensions defined in Section 2 have the 
values: 

ag = dy =r9(2/3)}, 
bo = (4n+ M)/X ](70/3?), 
co=[2(n—m)/> ](r0/3?), 
eéo=((4m+M)/> ](70/3), 
sin a= (2/3)?; 
sin Bo= 23/3; 


(30) 


COS a= 37; 


cos Byo= — 1/3, 





where >> denotes the total mass, }> =2m-+2n 
+ M, of the molecule if m is the mass of H, » the 
mass of D, and M the mass of C. The kinetic 
energy coefficients defined in Eq. (4a) become: 


My =(2m+M)(2m/M); po=(2n+M)(2n/M); 


we=(4um/M); p3=2m; py=2n; 
Ms=(m-+n)(2n/m); we=mitm; (31) 
Mz = Mot (n?/m) ; Me7 = Mi2— 1; 

Ms=mit(m?/n); wo=MetN; bso = mre — 7. 


The coefficients defined in Eqs. (23a), (24a), 
(25a), and (26) become: 


A,= —27343;=371(M+2m)/M: 
Ag= —2-3Ay= —3-}(2m/M); 
Ayp= —By=67}; 
B,= —27B;=3-1(2n/M); 
B,= —27°By= —3-"(M+2n/M: 
C= —D,= — (2/3); 
C;= —D,= —3-}; 
C,.=D,=C,=D3=0; 
Co= —Do= —373/2; 
A,=—3(m+n)/2m; 
(Ag we) = (Be, wer) = (2/3)!/2m; 
A;=34[1 —(ue/3m) ]; 
By= —31+(n/m) + (2yer/3m) 1/2; 
(As/uso) = (Bs/ uo) = — (2/3)? /2n; 
Ag= 3411+ (m/n) + (2ys9/3n) 1/2: 
By= —3![1 —(yup/3n)]. 


(32) 


For CHD, all the valence force constants ¢; 
carrying a given subscript in Eq. (28) become 
equal and the general valence potential function 
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Up is given by: 


2Uo=c1 T(6rm)?+e2r6? E (58a)? 


m ™m, % 


+63 i 67,67 m+ C4%o ps 5B mn (OF m +6r,) 


m>n m,n 


+Cs¥o i 5B mn(O% m? +47 n) 
+ cere : 6B mnOB m’n’ 
+¢77%¢" 2. re" (33) 


Relations between the valence force constants 
of Eq. (33) and the K’s of Eq. (27) are obtained 
from Eq. (29) simply by removing the primes 
on the c’s. 

It can be shown that the generalized force 
constants of Eq. (5) are the following functions 
of the valence force constants, for the case of 
CH.D,: 
ky _ [A + (2u1 ‘M)B), 3, 
keo=[A+(2po/W)B)/3, 
ky»p= [C+2(M+n4 wy) (B, M) | Ss, 
k3=kg=[40, + 500+ 2034 23(c4 —c5) 

+ (¢5/2) —3e;], 3, 
(34) 
ks = [ — Cy — 403+ 23(c4 - Cs) ba (5¢6, 2)+3e; ], 5 
ky=2'7D+(2m/ M)E)/3, 
ky=20F+(2m/MDE)/3, 
ko = 2'*7F+(2n/M)E)/3, 
koa = 237 D+(2n/ M)E)/3, 
where 
A = 2¢,+1002+ ¢3 — 2 -23(c,— 5) +66 — 6c, 
B=24+8c2—c3—4-23(cs—cs5) — 46, 
C= —2¢.—2¢3—2-23(e4—c¢5) —5e6+ 6c, 
D= 26, — 502 +¢3— 2 (cy — Cs) a (C6 2) + 3c;, 
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E = 2¢, — 40, — 63 — 27364 — 5) +266, 


P= o~ 2¢3 i 2-*e, aie Cs) + (5G¢ 2) — 3c: 





and 
kg =6[ (m+n) /2m P(2c2+¢6—2c), (35) 
kgp=Ae?(2e, —€3) +2A7*(2C2— Ce) 
+4A¢A7(C,—6s), 
k; = Be?(2¢e, —c3) +2Bz*(2c2— ce) 
+4B.Br(c.—Cs), 
ke; = AgBo( 2c, — 3) +2A7B;(2c2 — Ce) 
+2(A6B;+A7Be)(cs—6s), 
kg=As?(201 — C3) + 2A9?(2c2 — Ce) 
+4A sAo(ey—65), 
kg = B?(2¢; —c3) +2By?(2c2— Ce) 
+4BsBo(cs=cs), (37) 
Reg = AgBa( 2c, — 3) + 2A9Bo(2c2— Ce) 
+2(AsBy+AoBs)(c1—65), 


(36) 


where Ag, Az, ete., are quantities defined in 
Eq. (32). It is obvious from Eqs. (34)—(37) 
that the 17 generalized force constants of Eq. 
(5) can be expressed for the case of CH.D,» 
as functions of 6 independent valence force 
constants,” namely, €1, C2, C3, C1 —C5, C6, and Cz, 
which occur in the most general physical poten- 
tial function of the valence type. It follows that 
it should be possible to determine all six of the 
valence force constants from the infra-red and 
Raman spectra of CH2De since this molecule 
has nine different vibration frequencies of which 
all are Raman active and all except one are 
infra-red active. 


15 As a check on the expressions (34)-(37), the secular 
determinant (6) has been set up for the case in which 
m=n and it was found to yield correctly the expressions 
for the CH, frequencies. 
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The rotation-vibration Hamiltonian, complete to second order of approximation, is set up 
for the bent XYZ molecular model. The allowed energies are calculated and expressed in term- 
value form, E=hc(G+F); the vibrational term G is given explicitly and the elements of the 
secular determinant are given for evaluation of the rotational term F. The valence-force form 
of harmonic potential function is discussed for the bent X YY’ model and normal frequencies 


of HDO are calculated. 


I. INTRODUCTION 


HE rotation-vibration spectra of many 

polyatomic molecules, when observed with 
the excellent prism-grating spectrometers now 
available, exhibit fine structure and other features 
which cannot be adequately interpreted on the 
basis of the simple theory in which it is assumed 
that the vibrations of a molecule are simple 
harmonic and independent of the rotations. In 
order that the spectroscopist may obtain the 
maximum amount of information concerning 
molecular structure from the spectral data, it is 
necessary that an adequate theoretical study be 
made of the dependence of the rotation-vibration 
energies on molecular configuration, force con- 
stants, interactions between rotations and vibra- 
tions, etc. Such studies'-“" have already been 
made for several types of molecules and the 
general methods involved have been discussed by 
Dennison,’ Nielsen," Wilson and Howard," and 
Wilson." 
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Triatomic molecules, particularly those con- 
taining H or D atoms, afford a good field of study 
on account of the relative simplicity of the 
molecular configuration, small number of modes 
of vibration, possibility of resolving. the fine 
structure of the bands, etc. The bent symmetrical 
XY-2 or H,O type molecular model has been 
treated quite completely by Shaffer and Nielsen* 
(SN) and subsequently Nielsen'® has interpreted 
the fine structure of several HO bands. The 
linear symmetrical XY2 or CO. type molecule 
has been treated by Adel and Dennison? and the 
linear XYZ or HCN type by Nielsen." A com- 
plete study of the bent XYZ molecule, to the 
order of approximation required for analysis of 
the fine structure of infra-red spectra, has not 
been reported in the literature. Studies of the 
normal modes and force constants have been 
made by Salant and Rosenthal,!? Van Vleck and 
Cross,'* and others; Libby!® has studied the 
purely vibrational part of the problem including 
anharmonic terms. It has seemed advisable, 
therefore, to make for the bent XYZ model the 
same type of calculation as that made by SN? for 
the bent XY» model in order that the basis may 
be laid for adequate analysis of the spectral data 
available on such molecules as HDO, HDS, 
HDSe, and the like. In this paper the vibration- 
rotation energy expressions are obtained, to 
second order of approximation, for the general 


16H. H. Nielsen, Phys. Rev. 59, 565 (1941); 62, 422 
(1942). 

17E. O. Salant and J. E. Rosenthal, Phys. Rev. 42, 812 
(1932). 

18 J. H. Van Vleck and P. C. Cross, J. Chem. Phys. 1, 
357 (1933). 

19\V. F. Libby, J. Chem. Phys. 11, 101 (1943). 
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should be interchanged. 


bent XYZ molecular model of which HDO type 
molecules form a special case. 
II. EQUILIBRIUM CONFIGURATION AND 
COORDINATES 

The assumed equilibrium configuration of the 
general bent XYZ model, shown in Fig. 1, is a 
triangle of height Jo and base ko+lo. The X 
atom of mass m, is situated at the apex of the 
triangle, and the Y and Z atoms of masses mp» 
and m3, respectively, are situated at the base 
corners. In order to specify the equilibrium 
(vibrationless) positions of the atomic nuclei, a 
preliminary right-handed body-fixed rectangular 
coordinate system éf is adopted with the 
molecule lying in the & plane, the coordinate 
origin at the center of mass of the molecule and 
the ~ axis parallel to the base of the triangle as 
shown in Fig. 1. The equilibrium coordinates of 
the various atoms are shown in Table I, where 
M=m,+me2+m3. 

A set of principal axes of inertia xyz is adopted 
for the equilibrium configuration in order to 
facilitate the investigation of interactions be- 
tween rotation and vibration. The xyz coordinate 
system, shown in Fig. 2, has its origin at the 
center of mass, its z axis coincident with the ¢ 
axis, and its xy axes in the plane of the molecule 
and making an angle e with the & axes. The 
equilibrium x and y coordinates are given by 


Xo0i:= Eo; Cos e+ 70; Sin €, 


(1) 


Vai=— £0; sin €+ 0; COS €, 


where € is so determined that the equilibrium 
value of the product of inertia, 


Lecy= Doi MiXoiVoi, 


vanishes (I.y:=I.22=0, because 20;={o;=0). It 
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Fic. 2 
follows that 
tan 2e=(A/B), (2) 
where 
A = 2m ho(ms3ko— malo), (2a) 


B=m3(m,+mae)ko? +me2(m+ms)lo? 
+2momskolo—my(motms)ho?. (2b) 


The equilibrium values of the principal moments 
of inertia are 


Te2=(Uez-T1')/2, 
Tey=(Lee+T’)/2, (3) 
T,.=[B+2m,(mo+ms)ho? |/M 

where J’= +(A?+B?*)!/ A. 


The principal moments of inertia of the HDO 
type molecule can be obtained from the above 
general expressions if one sets ko=/p. 


III. THE NORMAL COORDINATES 


The instantaneous positions of the respective 
atoms during vibration are referred to a body- 
fixed right-handed rectangular coordinate sys- 
tem xyz, whose origin remains at the center of 
gravity of the molecule, whose xy plane is the 
plane of the molecule, and which coincides with 
the principal axes when the atoms are in their 
equilibrium positions. The instantaneous position 
of the zth atom is given by (xj, yj, 2:;) = (xoi +x’, 
Yoit+ yi, Zitz), where the primed quantities 
are rectangular components of displacement from 
the equilibrium position and where 20;=2,’ =0. 
One can readily choose a set of 3N—6=3 
generalized coordinates, suitable for the study 
of both the purely vibrational problem and the 
interactions between rotation and vibration, with 
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the aid of the Eckart conditions :'**" 


Lim«x'=0, Yimiy/=0, 


i m(xoiVi’ — yo,’ —Voix;’) =(). 


and 


An appropriate set of generalized coordinates, 
which reduce in the limiting case to those em- 
ployed in the XY:2 problem,* !7*! is 


Uy =X1' —[(meoxe' +m3x3')/(me+ms) |, 
U2 = yi’ —L(maye’ +msgys') /(me+ms) ], 


U3=Xo' —X3'. 


(4) 


It follows from the Eckart conditions and Eqs. 
(4) that 


xy’ =(p/m,)u4, 

x2! = —(u/(m2+ms3))ui+(u'/me)us, 

x3’ = —(u/(mo+ms))uy— (yu /ms) us, 

yi’ = (n/m) Uo, 

yo’ = (wa/mes) uot (u/s) +(u'y M2) Us, - 
ys’ = (uB/ms3)s— (pu! /ms)u, — (p’y/ m3) Us, 

















TABLE I, 
mt £0i Nor oi 
X m 1° (msko—moalo)/M (mo+ms)ho/M 0 
Y mz 2 (mskot+(mi+ms)lo)/M —mho/ M 0 
Z ms 3 (—moelot(mitme)ko)/M —myho/M 0 
where 


u=(m2+m3)m,/M, 

hp’ =m2m3/(m2+ms), 

ue” =myyor/ (X02 —Xos), 
a = (X03 — X01) / (X02 — X03), 
B= (X02 —X01) / (X02 —X 03), 
¥ = (Yo2—¥os)/ (x02 — X03). 


The kinetic energy of vibration of the XYZ 
molecule is 
T= (1/2) [msrtes? + pootte? + wsgtts? + 2urottitde 
+ 2pristtrtist+ 2postiotts |, 


mo/or 


Min=et(u’?/y’), 


(6 

where ) 
Moo = (u?/my) + (ua? /me) — (u?2B?/ms), 
uss=u'(1+7?"), 

20C. Eckart, Phys. Rev. 47, 552 (1935). 


* W. H. Shaffer and R. R. Newton,-J. Chem. Phys. 10, 
405 (1942). 


AND 
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(wi2/m’’) = (mos wy) = (ua me) —(uBms), 
Mis=ph’’y. 


The harmonic (quadratic) potential energy 
function is 


Uo=(1 2) (Riis? +Rootto? +333” 


+ 2R ott Ho +2kh 130 3+ 2ko3t ots |, (7) 
where the coefficients kj;, etc., are generalized 
force constants; the special valence form of har- 
monic potential function is discussed in Section 
IV for the case of HDO. 

The characteristic frequencies w, of the normal 
modes of oscillation are given in wave numbers 
by 


w= (24c)—dr,?; k=1, 2, 3, (8) 


where the , are roots of Lagrange’s secular 
determinant, 





| 

MiA—Rir | MidA—Rie | wish—kes 
MirAvA—Rie | boedrt— Ro» | MosA— Ros |=0. (9) 
Mist —Ri3 | Mo3d— Ros | Mast —R33 











The intermediate coordinates, “1, u2, and uz, 
are given in terms of the normal coordinates Q, 
by the relations, 


Uy =N101 +2202 +11303, 


U2 = N10 +N2202-+ 2303, (10) 
U3 = 23101 +N3202+N3303, 
in which the coefficients 1; are 
Nj. = Njx/ Ne, (11) 


if Nj, is the cofactor of the jk element in 
Lagrange’s determinant (9) with A=d,, (k=1, 
2, 3), and 


1 
2 


3 
m={¥ aN +E bi NaN || « & fj. 
j=l i’ 
The coefficients nj, can be determined for a par- 
ticular molecule from the experimental values of 
Ax, if there is reason for assuming a potential 
function Up») which involves only three force 
constants or if six force constants can be evalu- 
ated from additional data supplied by an isotopic 












































SPECTRA OF BENT XYZ MOLECULES 507 


molecule; the exact form of the normal modes The energy associated with the normal modes of 

can be found by combination of Eqs. (5) and vibration has the Hamiltonian form 

(10), in a manner similar to that?! employed for 3 

H.O and DO. IIo, = (he/2) ¥ wx (p.2/h*)+q.2), (12) 
Dimensionless coordinates g,. defined by the _ 

relation q.=(2mcw,/h)'Q,; are convenient for where p,=07/dq, is the momentum conjugate 

the quantum-mechanical part of the problem. to q;. 





The cubic and quartic portions of the potential function contain the following terms contributing 
to the first- and second-order energies: 
U; = hc[B11191° +8 22292* +833393° +811291°G2 + B11391°93 + 822192791 + 82239293 
+ 833193°i1 t+833293°d2+6123919293), (13) 
Uz= hcLyu1q '+-y 900092! + 333393! + 7112291792" + Y113391°93" + Y 223392°93" |: (14) 
IV. VALENCE POTENTIAL FUNCTION FOR HDO 


The HDO-type molecule is the special case of the bent XYZ model in which X becomes O, 
Y becomes H, and Z becomes D. It is reasonable to assume that HDO has the same equilibrium 
geometrical configuration as H.O and D.O, viz., an isosceles triangle with O at the apex. If Bo =2ap 
and ro denote, respectively, the equilibrium values of the apex angle and OH (or OD) distance, it 
follows that ko=lop=po/2=rosin ao and ho=rocosay. If r2=ro+dre, 73=7o+4r3, and p=pot+sp 
denote, respectively, the instantaneous values of the OH, OD, and HD distances, appropriate dis- 
placement coordinates for setting up a harmonic potential function of the central-force type are 
dro, 5r3, and 6p. If B=B +68 denotes the instantaneous value of the apex angle, appropriate coor- 
dinates for a valence-force potential function are 672, 673, and 68. The displacement coordinates 
defined above are, for infinitesimal vibrations, the following functions of the coordinates 1, #2, and us3: 


672 =A ju, +A ote +A gus, 

6r3= Byu, + Bot2+ Buu, 

68 = (Cy + Cous+Czus; }/ro, 

6p =D ym. +- Dou2t+ Dus, 
A ,=([(x01—Xo2) — (u’”/me)(Yo1— Yor) J/ ro, 
Bi= [(xo1 —Xo3) + (u’/ms)(¥o1— Yos) J/To, 
D,=(u'"/u’) (Yor — Yos)/ Po, 
Ci = (ro/ho) D1 — (po/2ho)(A1+ Bi), 
Az=p(yor — Yoo) (1/my) = (a/mz2) }/ro, 
B= p(yoi— oa) (1/1) +(8/ms) J/ro, 
D2= p(Yo2— os) (a/me) +(B/ms) ]/ po, 
C2= (r0/ho)D2— (po/2ho)(A2+Bz), 
A3= —(u'/me2)[ (x01 — X02) +Y(Yo1—Yo2) J/ To, 
Bs=(u'/ms)( (x01 —03) + (¥o1 — Yo) 1/70, 
D3=[(x02—x03) +¥(Yo2— os) |/Po, 
C3= (r0/ho)D3— (po/2ho)(As+Bs). 


where 


The complete valence-force potential function Uy for HDO is 


Uo= (1/2) [e1((6r2)? + (673)”) +c27 0" COs *axo(5B)? + 2cgr9 COS a5B(5r2 +6573) + 264572675 ]. 
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The generalized force constants in Eq. (7) are given in terms of the constants in Eq. (16) by the 
following relations, in which the indices k and k’ independently assume the values 1, 2, and 3: 


Rin = | (A;Ay +BpBr at Cy. Cr C2 COS 2a +[(A r+Bx) Cir + (Ay +By) Ci. les COS ao 
+(A,By +A, B,)cq}. (17) 


The constants in the valence potential function and the normal frequencies of HDO can readily 
be calculated with the aid of Eqs. (9) and (17) if one uses the force constants computed by Shaffer 
and Newton”! and molecular dimensions found by Darling and Dennison” from the spectra of H2O 
and D.O. The force constants, in units of 10° dynes per cm, are c; = 8.423, co= 2.049, c3;=0.4047 and 
c,= —0.1015. The molecular dimensions are 77>=0.9580A and By)=2a9=104° 31’. It follows that 
the angle e defined in Eq. (1) has the value e= 21° 5.5’ for HDO, and xo, =0.07093A, xo2=0.5667A, 
Xo3 = —0.8469A, yo1=0.07251A, yoo= —0.7472A, yos= —0.2020A. The kinetic energy constants 
defined in Eq. (6) have the following values for HDO (H =1.00813, D=2.01473, O=16.00000) in 
units of 10-*4 gram: w11= 5.886, woo= 5.814, w33= 1.282, wie= — 1.628, wi3= —0.5256, and w23=0.5142. 
The potential energy constants defined in Eq. (7) have the following values for HDO, in units of 105 
dynes per cm: ky, = 22.67, koo= 20.56, k33= 3.934, kig= —13.86, kis= —3.652, and ke3= 6.873. Sub- 
stitution of the above values of constants into Eq. (9) leads to the following values of \,, in units 
of 10” sec.—!: Ay = 5.3470, Ao = 2.8175, \3=0.77551. The calculated values of the normal frequencies 
of HDO, w,.=d,?/2mc, are w;=3882.(2) cm—!, we=2818.(1) cm, and w3=1449.(6) cm='. Libby,’ 
using slightly different values of the valence force constants given by Darling and Dennison,” has 
computed w;= 3883.8 cm™'!, w2= 2820.3 cm-!, and w3= 1449.4 cm~!. Two observed infra-red bands 
assigned”* to HDO have centers at 2720 cm and 1403 cm™. 

The constants computed above enable one to compute the following values of the principal 
equilibrium moments of inertia of HDO, in units of 10-*° g-cm?: J,,=1.210, I.y=3.070, I..=4.280. 


V. THE CLASSICAL ROTATION PROBLEM 


The principal axes of inertia for the equilibrium configuration of the bent XYZ model were 
specified by Eqs. (1) and (2), and the equilibrium values /,,, Z.,, and J,:, of the principal moments 
of inertia were given in Eq. (3). The instantaneous values, during vibration, of the moments and 
products of inertia are: 


3 3 7 
I= Ta 14E a Qetange+ LC ake gg ; 


k<k’=2 J 


k=1b k<k’=2 i 
(18) 
3 ? 3 , 

I,= In| 1+ 7 Cx’ Od +n gir+ 7 Crea ied’ 
k=1b k<ht=2 J 


ir 3 7 
y= In +E bi! det binge + ik bi Qed’ |. 
| 


5 








f j 3 
a | ’ , 
Izy= Teel} ) » dy’ Qe dingn? + DL dex Qe Qn 


k= k<k’=2 4 


——— 
~~ 


_ 


in which, if r=(h/2zc), 
a = 04a,’ =47-1A eL 8m Wotu+te(Yortayos Byos) + 2u’¥(¥o1 = Vos)” 3h }. 
by = wibby’ =47 3B. m xo tu’ (X02 — X03) Max], 
Ce = nex’ = C.L(ai/Ae) + (0;./B-.) ], 
dy =ox'd,’ =47 1A SBS Lmiyoinirtu’ (yoo — os) M3x |; 


2B. T. Darling and D. M. Dennison, Phys. Rev. 57, 128 (1940). 
*3Ta-You Wu, Vibrational Spectra and Structure of Polyatomic Molecules (Prentice-Hall, Inc., New York, 1940). 
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» _# , 9, 
Ak = CK POI? A jek? > 2A .[ (u"’ / u’)n UM 1k HH M22QN aN oK" te MEE 


+ o12( MyM + 1M 1%) + p23(Mox.N3x° +N yMNoK") + is( MLN sK + 13.N 1%") |, 


1 1 U 
Dye = 2K? Dj = 2B, [un Mike tu’ NgiN3x |, 


onlay ~2C + Sie = CL (aus A)+(bix/B.) ], 


kh’ ¥k 


Chk 


Cone = neon 8px F4C, Zz. Carer Serer = CL (arn A d+ (dix /B.) |, 
Re? 


dine =2A PBA ( my ne + neni) Fe” (MmyMgice + get?) +e” (MoxMgne + Ngee) + 2p’ yng Nzi-_|, 
Cate = — Cher = MM yMogr — NM) Fe” (MgMee — NrgNgn) Fy” (MgeMoxe — NxM 3x"), 

=m 1Xo1/(Xos — X02), 

and the equilibrium values of the rotational constants are 


A,=(h/8r'cl,,), Be=(h/8r'cl,.,), C.=(h/8rcl,2). 


The quantities a,,, b,;, and d,, are obtained, respectively, from a,x, bx, and dy. by setting k’=k. 
It follows from the definitions of c; and c;,, that 


(¢./C.) —(ar/Ac) —(b./B.)=0 and (Cyr / Ce) — (dine /Ac) — (dix: / Be) =0. (20) 


The vibrational angular momentum of the bent XYZ model has only a z component, different from 
zero, 


P2= Do Seer L(wer/ wr) *QuPur — (we/ wer) *Qurpr J. (21) 


If the rotating vibrating XYZ molecular model is considered as a semirigid rotator with an 
internal flywheel, its rotational energy is given in the classical Hamiltonian form by 


Trot = (1 2) {Mere tH hyyPy? + uee(P:— p:)?+2uryPP,y}, (22) 


where P,, P,, and P, are components of total angular momentum, p, is the internal angular mo- 
mentum, and, to second order of approximation, 


i , ° 9 
Mer=Ler {1 — Qo Lan’ ge — (ae? — an tds?) qi? }}, 


k 
by = Tey (1 —X [bulge — (be? — Bie de?) qu? J}, 
k 


1 » ’ 9 
Mez=Lez {1-0 Len’ ge — (x? — Cin) gu? J}, " 


k 


$4 
Mery = "a | > (dx qs + (d’, k —a,'d;’ _ by’ dx’) qi? | } . 
k 


(Quadratic cross-product terms have been omitted from the above expressions because they do 
not ordinarily contribute to the first- or second-order energies.) 


VI. THE QUANTUM-MECHANICAL HAMILTONIAN 


The appropriate quantum-mechanical Hamiltonian for investigation of the rotation-vibration 
problem of the bent XYZ molecular model can be set up according to the method of Wilson and 
Howard" and has the following zero-, first-, and second-order parts which can contribute to the 
allowed energy expressions to second order of approximation : 


3 
Io = (he/2) D wl (pi2/h?) +442] + (he/h2)(A P2+B.P,2+C.P."), (23) 


k=l 
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3 
Hy =U, —(p:P:/Iez:) —(he/h?) © {alan AcP2+b,’B.P,y?+c’C.P2—d;'APtBH(P:Py+P,P:)}}, (24) 


k= 
: 3 


H2= U2+(he/h®)[Cep2)+(he/h®) E bq? (an’?— arn tdy!2)A P22 + (by? — Bin +d) BP? 
k=l 


+ (cx? — cin) C.P22+ (dix, —ax'de! —by!dy!)A EB (PP, +P,P.)]}. (25) 


Certain lengthy terms, which usually occur in JJ; and H2, have been omitted from Eqs. (24) and 
(25) because they either have no effect on the energies to second order of approximation or affect 
only the constant term /cGp in the final energy expression (28). 

The perturbation calculation required for the evaluation of the eigenvalues of the Hamiltonian 
is simplified by a contact transformation, H’=T7THT-', of the type described by Thomas** and 
others,*"°" where T=e®S (i=4/—1) and Hy’=Hy, Wy’ =H,+i1(SH»—H)S), and H,’=H, 
+ (i/2)[S(i1+H1’') —(Hi+H1’')S]. In the cases where there are no vibrational degeneracies, real 
or accidental, it is possible to set up a function S which reduces H,’ to zero and the perturbation 
calculation involves only the evaluation of the eigenvalues of Ho’ and H,’. For the bent XYZ model, 
H,! has the following form (each of the subscripts k, k’, and k’’ can assume the values 1, 2, and 3, 
but in a given coefficient kk’, k'*k’’, kk): 


Hp! = hcl p> | veo +E vhawvaran? +E (qx?/h?) [aA -P.?+8,B.Py+y.CP 
J ke i 


k 


+5,A3B.(P,Py+P,P:) ]+(1/h2)(C.p2)} +H (26) 


where H2,, the second-order energy correction associated with centrifugal stretching, has the form 
Hoe= = (he/h*)[trr22P 23+ faak + Tessol ot Teal Pele t+ PPA) + teed PPA +P fP,*) 
+ tyye2e( PP 2+ P2P,?) + teyzy(P2Pyt+P,P:)?— (tend o' + Tunet o° 


f + Te2ryP2*)(P2:Pyt+P,P:) —(P2Pyt+PyP)(t22r2yP 2+ Tyyysl yy oe a BF (27) 
1 


YikKh = Ez _ (5Binx/2cx) nae > (8c, hans 3u?) (Benn: / eon) / (Aeon? ees wy?) | 
ke 


veenrer = [Yann — (QeonBirx-/ (4e,.? — wpr?)) — (2x Bernre/ (Awe? — 4?) — (3B ci Benn /@x) 
— (3B xn Burnrn/wnr) — (Berne Buren / wn) — (wo — wy” — wy”) (cox 128/2D125) ]; 
Dy23= (@1+w2+w3) (w1 — we — w3) (W1 — w2-+ 3) (wi to2— ws), 
ore = (4? — anny”) /ewr) + (SaxBene/ wx?) + . (AxBrrn/wx?) |, 


B= ( (du? — dix +d”) /wox) + (3b: Been / wx!) + a (de Bren? /wx?) 1, 


ve =L((cx? — Cun) /oe) + (36:Bire/OB)AY (CeBrrnr/wx®)+2C. Y (Sap? + wx?) OE eas — 
bx = ( (dus, — and, — byd;) /w. }, k | : 
Trrez=Ac® Qin (ax?/2ux”), Tyyw = Be? Lx (be*/2ax”), 


Te222= Ce? Dix (C4?/2wx?”), Trzyy=AcBe Dx (axb;/ 20”), 
Trrzz=AcCe Dix (axcx/2wx”), Tyyz2e=BeCe Dex (diCx/ 20%”), 
Tayzy=AcBe Dx (dy? /2u,”), Trzxey= APB Dix (andy /2u”), 
Tyy2=APB) Dix (didi /2u7), Tezexy=ABACe Dix (Cedi /2wx”). 


Inspection of Eq. (26) shows that resonance degeneracy can occur in any of the following cases: 
wp wr, wwe tw, OF ww. In such a case, the term from which it arises cannot be removed 


41. H: Thomas, J. Chem. Phys. 10, 532, 538 (1942). 
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(24) from HH,’ by the contact transformation and a degenerate perturbation calculation is required to 
evaluate its effect on the energy. 


VII. THE ALLOWED ENERGIES 





(25) The allowed values of Hy and Hz have been determined by a method similar to that employed 

and | by SN? for the bent XY» model. The first step in the procedure is the evaluation of the general matrix 

fect element of Ho’+H.’ which is diagonal in all the vibrational quantum numbers, while the angular 
momentum operators are considered as constants. The result can be written in the form Eyi,+/7e, 

Man where Eyi,=hcG contains all the terms in the energy which are independent of rotational quantum 

and | numbers and Hz is a form of the Hamiltonian for an asymmetric rotator with the equilibrium rota- 

=H, tional constants replaced by values effective in the respective vibrational states and with the addition 

real of centrifugal stretching terms. 

tion The vibrational term G is given by 

- G=GotL Con(v. +3) +Gru(ve +3)? + » Gixe(Ve +3) (de +4) J (28) 

where , _ 


Go= (1/16) Ss [Gyasse— (7Bane/con) + (Seon-Beae’ / (4e0n® — con®) — (4eor0020008122/D 192) ] 


+ (1 /4) 7 | au -tbu-ben-babibhierben —d,?— (3/8) (a. +b. +x)? —2C. = | 
k 





(26) z= 
_ Gus (1/4) [6ys111~ 15 Bins/«1) ~ Bxa(Bin2/e2) ~ Exa(Bins/es)} 
Ga2= (1/4) [672222 — 15 (B222/w2) — E21(B221/w1) — E2s(B223/0s) ], 
| Gaa= (1/4) [673333 — 15(8333/00) - Exi(B3s1/1) - Es0(Bsa2/w2) ], 
27) | Gy2= 1122 — 3(B1118221/w1) — 3(B2228312/w2) — (81138 223/ws) 


. — F12(B512/1) - Fox(B201/«2) - Cs12(Bi28/ws) 4+Kish isCor 
G23= 72233 — 3(B2228'332/we) — 3(83338223/w3) — (Bs31B221/w1) 
— F3(B 223/02) — Fs2(Bss2/10s) —Ci29(Bi23/01) + Kash 29Ce, 
3) J, Gai = 3311 — 3(83338331/3) — 3(81118113/1) — (81128 332/w2) 
— F13(B 313/01) — Fu(B3s1/ws) —Cos1(6128/w2) + KatnCe, 
Eur = ( (8c? — 3c?) / (400? — wn?) J, 
Fyxs = [2042 / (40%? — wx?) J, 
Creer = Lon? (wx? — oye? —wy/2D123)), 
Kx = [(wx? +e?) /wxorns J. 


The effective rotational Hamiltonian Hp is 





if 


Hp = (he/h*) {AyP.2+ByPy?+CyP2+Dy(P:Py+PyP2)} +H, (29) 
; where H;, is the term defined in Eq. (27) and 
Avy=A{1+ Xi on(u.+3)], 
By =B(1+ Lx Be(u+3)], (29a) 
es: by Cy=C.L1+ Xk ve(ue+3) 1], 


ved 4 


Dy=A3B} ¥ x 6(vx +3), 


if ax, Bey Ye» Sky Terre, etC., are quantities in connection with Eqs. (26) and (27). 
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It follows from Eqs. (29a) and (20) that the quantity A=[(1/Cy)—(1/Av)—(1/By) ] is inde- 
pendent of the anharmonic constants, to second order of approximation, and has the value 


A — ar ween? Lee -, (b,?+d;”) I, — (a,?+d,.”) I, m3 (30) 


A similar relation has been found to hold for the bent XY» model by Darling and Dennison and 
for the planar XY; model by Shaffer and Silver ;° such a relation does not appear to hold for non-planar 
molecules.® 

The secular determinant for evaluation of the allowed values of J7,p is similar to the Wang” secular 
determinant for the asymmetric rigid rotator, and like that proposed by SN* for the semirigid 
bent XY» model. If the wave functions for the symmetric rigid rotator are used as basic functions 
for the expansion of the rotational wave functions of the bent XYZ model, the secular determinant is 


|(JKM|Hp|JK'’M)—Erdx| =0, (31) 


where J, K, and M are the usual quantum numbers associated, respectively, with total angular 
momentum [P?=J(J+1)h?], component of angular momentum along the rotator’s symmetry axis 
[P.=Kh], and component along the direction of a magnetic field [Pr = Mh]; J is a positive integer, 
K=0, +1, ---, +J, and M=0, +1, +2, ---, +J. The appropriate matrix elements for setting 
up the secular determinant (31) have been given by SN*, are all diagonal in J and M, and have 
values independent of M. It is necessary to solve a 2/+1-fold determinant for each value of J with 
the elements, 

(K | Hp | K) —Er=he { Rit R2K7P+R3K4} _ Er, (31a) 
(K |Hp| K+%2)=he|(RstiRy’)+2(Rs+iRs’)(K°?+2K42)}[f-—A(K41) JL f—(K+1)(K+2) }', (31b) 
(K|Hr|K+4) =he(RetiRe’)[f—A(K +1) PL f—(K+41)(K +2) }'Lf—(K+42)(K+3) }} 

x f-(K+3)(K+4) ]}',  (31c) 
where , 
i=(-1)', f=J(J+1), 
Ri= (f/2)(A v+By) +> (f/4)(tTa2r2+ Tyyuy - 2 tesyy FAT eyey ] oe (f?, 8)[3ter22+3 Tyyyyt 2T cyt 4 Teysy); 
R= (1/2)(2Cy —By—A v) = (5/8) t2222+ Tyyuy 2 recy t 4 Teysy | 
+ (f/4) [3 rerez t3 TyyyyA2 Terry —4 T2222 —4Tyy2z2 FA Tryzy |; 

R3= (1/8) [ —3recrz— 3 Tyyyy— 87 2222 — 2Tr2yy A B8T2222zF BT yy22—4Tryzry_]s 
R= (1/4)(Av — Bry) —(f/4)(tr2r2— Tyyun)s 
R= (1/2) Dy +f( Trezyt Tyeys) |, 
R;s= (1/8) rrex2— Tuyyy a | a 
Rs = (1/4)[27222y— Trxery Tyeve Js 
Re= (1/16)[— Vases Tyyyy t 27 csyy tA Teyzy_); 
Re = (1/4) [t222y— Ty: )- 


The complex character of the off-diagonal elements defined in Eqs. (31b) and (31c) is associated 
with the complete lack of point symmetry in the configuration of the bent XYZ model. The secular 
determinant for each value of J can be factored into four parts, in general, and expanded in the 
manner discussed by SN* and Nielsen for the bent XY» model. 
The total rotation-vibration energy of the bent XYZ model can be written in term-value form, to 
second order of approximation, as 
E=hc(G+F), (32) 


%S, C. Wang, Phys. Rev. 34, 243 (1929). 
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where G= Ey/hc is the vibrational term given in Eq. (28) and F= Er/hc, if Er is a root of the secular 
determinant (31). The selection rules allow all possible vibrational transitions because of the asym- 
metry of the molecular configuration. The selection rules for the rotational transitions and the 
calculations of the intensities of the rotational lines in the bands of asymmetric molecules have 
been discussed by Dennison,!! Nielsen,'® and others. 
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Electrical Anisotropy of Xerogels of Hydrophile Colloids. Part II*+ 


S. E. SHEPPARD AND P. T. NEWSOME 
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(Received October 5, 1944) 


An improved, though still relative, quantitative expression has been derived for the electrical 
anisotropy E. A. The E. A. is a linear function of the elongation up to a certain limit, as is also 
the optical birefringence. The relation of electrical to optical anisotropy has been studied in 
greater detail. A more critical discussion in terms of atomic model structures is given of the 
hypothesis that E. A. in the compounds studied is due to the formation of continuous par- 
allel chains of hydrogen bridges, having electronically conducting character. The materials 
studied were: polyvinyl acetate and its hydrolyzed stages down to polyvinyl alcohol, cellulose 


acetate and its hydrolyzed stages down to (hydrate) cellulose. 


N Part I the electrical anisotropy was ex- 

pressed quantitatively in a quite relative 
form as ¢/V7h, where ¢ is the angle of rotation 
(from 45° as rest) by a field V’, for a cylinder of 
unit radius and thickness /.** 


MODIFICATION OF THE QUANTITATI 

EXPRESSION 
It has been observed that the effect, for a disk 
not too large compared with the electrode, or 
pole-piece diameters, is an inverse linear function 
of the square of the distance between them, con- 
firming the proportionality to V*. Further, it is 
found that the effect, other things being equal, 
is proportional to r* where r is the radius (cf. 
Fig. 1). Thus, it is proportional to the volume hr’, 
of the cylinder or disk, and again to the radius 
r, as defining the torque on the dipole of effective 
charge separation 2r. To compare substances of 
different (monomeric) molecular weight, we can 

advance to the expression: 


E. A.=K-f(Al/DhrV7d/ Mn, 


where K is an apparatus constant which can also 





Presented in the inaugural program of the Division 


of High-Polymer Physics, American Physical Society, 
Rochester, New York, June 23-24, 1944. For Part I, see 
}; Chem. Phys. 12, 244 (1944). 

t+ Communication No. 997 from the Kodak Research 
Laboratories. 
** Instead of ¢, used in Part I. 
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Fic. 1. Variation of electrical anisotropy of polyvinyl 
alcohol with the radius of the disk (200 percent elonga- 
tion). 


include Avogadro’s constant, r=radius in cm, 
h=thickness, f(A//l) is a function of the elonga- 
tion, which reduces to proportionality for a 
certain range, V=voltage, d=specific gravity, 
and M,,=“‘equivalent’’ monomeric molecular 
weight.* 

~ *E.g., Cx(OH) for polyvinyl alcohol, Co(OH); for 


cellulose. 






































































































514 S. E. SHEPPARD AND P. T. NEWSOME 
1000 / 2000 | 
150 / 
1s00 
' 
2 
*e 
= 500 eo 
S y © L000 
| 
w 
ons /| ] 
500 
1 
, Vv 
° 100 200 300 iat] 
*/e ELONGATION 
, oes 4 


Fic. 2. Effect of elongation on electrical orientation of 
polyvinyl alcohol. 1. High viscosity. 2. Low viscosity. 


ELONGATION 


Since, usually, internal orientation or ‘‘fiber- 
ing” has been produced by stretching, the rela- 
tion of E. A. to the amount of elongation to which 
the material has been subjected is important. 
This will be discussed in greater detail later; at 
this point it suffices to state that over a certain 
range, and with certain conditions maintained 
constant, the E. A. appears proportional to the 
actual elongation Al/], which, usually, will be 
expressed as a percentage. This is illustrated in 
Fig. 2. Fuller confirmation of the proportionality, 
up to some kind of elastic limit, is shown by the 
relation of the optical birefringence to the E. A. 
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MOLES OF WATER PER MOLE OF HYDROLVZED VINYL ACETATE 


Fic. 3. Effect of water content on the electrical ani- 
sotropy of polyvinyl acetates at 225 percent elongation. 
1. Polyvinyl alcohol. 2. Hydrolyzed polyvinyl acetate (73 
percent polyvinyl alcohol). 3. Hydrolyzed polyvinyl acetate 
(35 percent polyvinyl alcohol). 4. Polyvinyl acetate. 





° to 20 
MOLES OF WATER PER MOLE OF CELLULOSE MATERIAL 


Fic. 4. Effect of water content on the electrical anisot- 
ropy of cellulose materials. 1. Cellulose, 0. 2. Cellulose 
acetate, 24.5 percent acetyl. 3. Cellulose acetate, 40.0 
percent acetyl. 4. Cellulose acetate, 44.5 percent acetyl. 


on the one hand, and to the elongation on the 
other (vide supra). These relations—of the linear 
regime to that beyond—are undoubtedly of 
great interest for the problems of structure and 
elasticity, but for the present we shall disregard 
them and confine ourselves to comparison of 
materials at corresponding stages of elongation 
—and in which, for the present, only “frozen in”’ 
states are considered. 


MATERIALS 


Previously, it was found that only colloids 
containing autonomic hydroxyl, or perhaps 
an equivalent amino-group, as in proteins, were 
reactive. It was further observed that polyviny! 
alcohol was reactive independently of moisture 
content, within measurable limits. At this stage, 
therefore, attention was largely concentrated on 
materials of the two macromolecular series: 


1. Polyvinyl acetate—polyvinyl] alcohol. 
2. Cellulose triacetate—cellulose. 


INFLUENCE OF ADSORBED WATER 
(MOISTURE REGAIN) 


Neither fully acetylated polyvinyl acetate, nor 
“cellulose triacetate,’’ shows more than a trace of 
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ELECTRICAL 


response at any humidity. Very slight residual 
effects at high humidity might be due to “‘spots”’ 
of incomplete esterification. With polyvinyl 
acetate there is practically no effect of moisture 
at either extreme of composition, but there is a 
definite sensitivity for intermediate, partly 
acetylated or deacetylated materials (cf. Fig. 3). 
With cellulose, while there is practically no effect 
of water upon the “‘triacetate,’’ there is a steadily 
increasing effect on deacetylating, no maximum 
being reached, but the greatest effect occurring 
with cellulose itself (cf. Fig. 4). 

Actually, the range of water adsorption—in 
terms of vapor pressure or relative humidity— 
which could be studied was limited by the 
rigidity of the sample. Above a certain value the 
disk becomes flabby and badly distorted. Even 
before this, the adsorption of water tends to a 
deformation, and there occurs swelling and dila- 
tion principally transverse to the direction of 
stretch, as shown in previous studies on gelatin.! 
The over-all behavior is therefore represented by 
that shown in Fig. 5 (for gelatin). The trans- 
verse expansion introduces a ‘‘form”’ anisotropy 
which first counteracts and then surpasses the 
internal anisotropy. 
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Fic. 5. Change of E. A. of gelatin with time at 95 percent 
R. H. due to unequal swelling along and across direction of 
elongation. Note: The open circle shown to the left of the 
above figure should be ¢. 


*S. E. Sheppard and J. G. McNally, Colloid Symposium 
cont a Wiley & Sons Inc., New York, 1930), Vol. 
,p. 17, 
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Fic. 6. Relation of E. A. to intrinsic viscosity (chain length) 
of polyvinyl alcohol samples. 


EFFECT OF MOLECULAR WEIGHT OR 
CHAIN LENGTH 


A number of polyvinyl alcohol samples having, 
according to measurements of their specific vis- 
cosities, considerably different average molecular 
weights, was examined. It appeared (cf. Fig. 6) 
that the E. A. is independent of the mean 
molecular weight and chain length, at any rate 
over a considerable range. 


CORRELATION OF ELECTRICAL ANISOTROPY 
WITH OPTICAL BIREFRINGENCE 


In Part I definite birefringence measurements 
were educed as evidence of “fibering,’’ of macro- 
molecular orientation, for the sheets examined. 
This relation has been studied in greater detail. 

As noted, with increasing elongation, the E. A. 
of polyvinyl alcohol increases regularly—and 
approximately linearly—up to about 250 percent, 
when the rate of increase of E. A. appears to 
decrease rather sharply (cf. Fig. 2). The optical 
birefringence was also observed to increase at 
first linearly with the elongation, but then a 
nonlinear region, with indications of approach 
to a limit, supervened (cf. Fig. 7). For the pri- 
mary linear relation between the E. A. and the 
optical birefringence for the same material, see 
Fig. 8. 

What has just been said applies only to the 
hydrophile colloids which manifest the property 
of electrical anisotropy, and not to the purely 
organophile ones, such as cellulose triacetate, or 
polyvinyl chloride which become birefringent on 
elongation, but not electrically anisotropic. 


Moreover, while in most cases the optical bire- 
fringence is positive, in the sense that n,>nz, 
where the maximum index n, is in the direction 
of the elongation, this is not necessarily the case. 
Examples were noted in Part I where n, could 
be perpendicular to the plane of the sheet, and 
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Fic. 7. Influence of elongation on the birefringence of 
polyvinyl alcohol and polyvinyl acetates. 1. Polyvinyl 
alcohol. 2. Hydrolyzed polyvinyl acetate (73 percent poly- 
vinyl alcohol). 3. Hydrolyzed polyvinyl! acetate (35 percent 
polyvinyl alcohol). 4. Polyvinyl acetate. 


it could happen that it would be in the plane of 
the sheet, but approaching a_ perpendicular 
direction to the axis of stretch. Such conditions 
occur, with higher alkyl esters of cellulose? and 
with polyvinyl acetate, when alignment of optical 
density, by side chains, may be greater transverse 
to the backbone spacings of the chain molecules. 
With partly hydrolyzed polyvinyl acetate, the 
optical birefringence ,—mg decreases with de- 
creasing —OH content from positive to negative, 
going through a zero value at about 6 percent 
hydroxy] (cf. Fig. 9). 


DISCUSSION 


The change from a well-marked linear relation 
(over a considerable range) to a non-linear one 
may well indicate something like a phase transi- 
tion, such as occurs in the “‘crystallization’”’ of 
rubber upon high elongation.’ In the case of 
polyvinyl alcohol, a sample at some 250 percent 
elongation gave such a well-marked fiber pattern 
(cf. Fig. 10) as to indicate that something of the 
sort is occurring, and this will be followed up by 
further x-ray diffraction observations. 

In the previous paper it was suggested that the 


phenomena observed might be accounted for by 
2J. Spence, J. Phys. Chem. 43, 865 (1939); 45, 401 
(1941). 
3H. Mark, Ind. Eng. Chem. 34, 1343 (1942). 


au Fr. 


NEWSOME 


the formation of oriented chains of hydrogen 
bonds, relatively rigidly connected with, or in- 
tegrated in, the long-chain structure of the colloid 
macromolecules. The outstanding effects which 
have to be explained are: 

1. Absence of effect with purely organophile 
(hydrophobe) colloids (polymers). It appears 
that Wood’s interpretation in terms of a differ- 
ential dielectric capacity along and across an 
axis of fibering, while a relevant, and, in a sense, 
obligatory concomitant, is not a sufficient cause. 

2. It is thought that the role of effective —OH 
groups, whether from constitutional 
hydroxyls (and perhaps amino- or imino-groups 
in proteins) or from these reinforced by adsorbed 
water molecules, is fo form quasi-electronically 


derived 


conducting chains. 

In the case of polyvinyl alcohol an idealized 
image of a parallel group of such chains, part of 
the molecules, is shown in Fig. 11. 

In the dry material it is supposed that the 
chain molecules are aligned, side by side, with 
alternate hydrophobe (—CHz:) and hydrophile 
—QOH adjacent alignment.t The extreme ease of 
lengthwise splitting supports this view. 
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Fic. 8. Relation between electrical anisotropy and optical 
birefringence of polyvinyl alcohol. 


A recent x-ray study of the structure of 
(crystalline) polyvinyl alcohol’ indicates that 
some modification of this picture may be neces- 


+ En bloc, as compared with a plane, some complications 
of the arrangement might occur. 
' Rose C. L. Mooney, J. Am. Chem. Soc. 63, 2828 (1941). 
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sary, at any rate for well-crystallized material. 
It was not stated how the material was prepared 
but merely that ‘‘well-oriented samples can be 
obtained, the fiber periodicity is short, and the 
chemical composition known and simple.”’ The 
structure concluded to represent the observa- 
tions best did not accord with the assumption 
just made that the hydrogens are lying (as 
bridges) between successive oxygens in the same 
plane. This would not explain the 2.72A separa- 
tion between oxygens of adjacent chains, and 
would give ‘‘a hydrogen bond distance somewhat 
shorter than is usual in alcohols. A more reason- 
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°/o HYDROXYL 


Fic. 9. Effect of hydroxyl content on the optical bire- 
fringence of hydrolyzed polyvinyl acetates at 225 percent 
elongation. : 


able arrangement is that illustrated in Fig. 2 
[our Fig. 12]. This shows paired vinyl alcohol 
chains projected in the 101 plane. The structure, 
as deduced from x-ray diffraction, is considered 
to consist of pairs of long-chain molecules, held 
together by hydrogen bonds. These double chains 
all lie along the fiber axis, but in random orien- 
tation about it. . . . Each chain has a zigzag 
carbon planar configuration, with glycol units in 
the 1, 3 position.” 

It appears that this paired structure would 
still permit of continuous hydrogen-bond chains. 








Fic. 10. Diffraction diagram, Cu Ka of polyviny] alcohol 
at 250 percent elongation, parallel to plane of sheet and 
perpendicular to direction of elongation. 


The ease of splitting longitudinally observed in 
our specimens would still be possible, provided 
the double chains were not in too random orien- 
tation about the fiber axis. Also it is possible, as 
already mentioned, that crystallization occurs 
rather abruptly as a function of elongation and 
orientation. Insofar as H-bonding and conduct- 
ance are concerned, the two molecular chains 
constitute only one shared H-bond chain. 

A kind of electronic conduction could take 
place along the macromolecular chain by way of 
a quantum molecular resonance process; the 
hydroxyl! can be regarded as resonating between 
the states symbolized below: 


—QH——O----H*t— —OH—etce. 


In an electric field a torque would be produced 
proportional to V?, to the total current capacity 
of the conductors, and to sin 6 cos @ («sin 26), 


\ : 

| F i ae H oH 
| | 

S's 7s 
Ny a — : : 
os Ff” ~- F F 
H o=H, H, = H H 
| 1 | 

7 va O-H 


Fic. 11. Formation of chains of hydrogen bonds in 
polyvinyl alcohol (adsorbed water may form a separate 
chain as shown). : 
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+ Fiber direction — 


Fic. 12. [From Rose C. L. Mooney, J. Am. Chem. Soc. 
63, 2828 (1941).] Polyvinyl alcohol chains projected on 
the 101 plane. The circles, in descending order of size, 
represent oxygen, carbon, and hydrogen. The dotted and 
solid circles distinguish between hydrogens on alternate 
sides of the chains. 


where @ is the angle between the axis of the 
conductor, and that of the electric field. The 
‘independence of length’’ for polyvinyl! alcohol, 
etc., seems somewhat at variance with the effect 
of adsorbed water with partly hydrolyzed poly- 
vinyl acetates. In this case, evidently gaps would 
exist in the molecular chains, consisting of 
residual acetyl groups. It appears that these can 
be ‘“‘bridged’”’ by adsorbed water molecules. 
Beyond a certain composition, the nearer this 
latter approaches 100 percent polyvinyl alcohol, 
the less effective are the additional water 
molecules. 

Considering again a completely stretched cel- 
lulose macromolecule or chain, then according to 
x-ray and other indications it appears that while 
the two secondary —OH groups in one glucosan 
unit can possibly supply an H bond between 
themselves, it would be quite difficult for a 
primary —OH of one glucosan to bridge to a 
secondary —OH of a contiguous glucosan. Only 
by adsorption of water molecules could the two 
chains of H bonds be formed, as compared with 
the one chain feasible for polyvinyl acetates (cf. 
Fig. 13). This seems a reasonable explanation 
why the E. A. of partly hydrolyzed polyvinyl 
acetates increases linearly at first with the 

5K. H. Meyer and H. Mark, Der Aufbau der hoch- 


polymeren organischen WNaturstoffe (Akad. Verlagsges. 
m.b.H., Leipzig, 1930). 


SHEPPARD AND P. T. 
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density of adsorbed water molecules, while with 
cellulose a power (initially 2) of the adsorption 
density seems involved. Exact relations are im- 
probable, particularly for cellulose, because 
thermal and other conditions would affect the 
strength of individual hydrogen bridges, or, con- 
versely, the interval necessary. Certainly for 
resonance and conductance along the chain, the 
strength of the H bond would have to be limited 
possibly up to 2.9A.° This could be satisfied 
readily enough in polyvinyl alcohol itself and 
probably without much difficulty in the ‘‘gaps’’ 
for hydrolyzed polyvinyl acetates, but in cel- 
lulose both the stereochemical and the energetic 





Fic. 13. Continuous hydrogen bonding in cellulose by 
the addition of one mole of water per glucose residue 
(bonds contributed by water are indicated by circles). 


conditions for the adsorption are considerably 
more complicated. The figure given (Fig. 13) is 
based on the configuration deduced in the paper 
of Meyer and Misch’ of 1937. In a recent paper 
by P. H. Hermans, J. de Booys, and C. J. Maan® 
certain modifications are suggested, which are of 
importance in regard to the regain of moisture 
(water adsorption) and to the bearing of this 
upon electrical anisotropy. These changes are 


®L. Pauling, The Nature of the Chemical Bond and the 
Structure of Molecules and Crystals (Cornell University 
Press, Ithaca, New York, 1940), second edition, p. 289. 
1937} H. Meyer and L. Misch, Helv. Chim. Acta 20, 232 

8 P. H. Hermans, J. de Booys, and C. J. Maan, Kolloid 
Zeits. 102, 169 (1943). 
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Fic. 14. [From P. H. Hermans, J. de Booys, and C. J. 
Maan, Kolloid Zeits. 102, 177 (1943).] (a) Stuart model 
of the cellulose chain with spatial puckering viewed from 
above (oxygen atoms shaded). (b) Diagram of the cor- 
responding nuclear structure. 


largely based upon a remodelling using van der 
Waals radii (according to Stuart) for the atoms, 
and instead of a planar pyranose (hexose) ring, 
a puckered one, according with the “armchair” 
form of the cyclohexane ring. 

In Fig. 14 is shown their ‘“‘most probable con- 
figuration of the straight, stretched cellulose 
molecule with spatial buckling, and _ using 
Stuart’s atom-models (van der Waals radii) 
viewed from above. It will be seen that the 
bridge oxygen between nuclei II and III lies 
somewhat downward, and, between | and II, or 
equally III and IV, somewhat upward [of a 
median plane] while the horizontal buckling 
factor is clearly evident. The model also shows 
that a horizontal buckling in the opposite sense 
would lead to steric hindrance.” * 

‘“‘Herewith only the choice of position of the 

CH:OH groups is left free. In Fig. 12 these 
[groups ] are placed alternately in the ‘normal’ 
and the ‘reverse’ positions. —OH group 3 is 
given a particular position with respect to the 
ring-oxygen atom of the adjacent ring, such as 
to give a hydrogen bond. Likewise, such an 
H-bond is possible between the —OH in normal 
position at C-atom 6, and the glucosidic bridge- 
oxygen. By these additional ring systems, the 
stretched position of the chain might be further 
strengthened, since rotation about the glucoside 
bridges would be hindered.” 

The authors also point out this notable peécu- 
liarity of the 8-glucose configuration [as adopted, 
cf. Figs. 2 and 3 of their article]. ‘All the hy- 


* Between hydrogens of successive glucosans. 
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droxyl groups find themselves on the side of the 
ring extension, while above and below only 
hydrogen-carrying valences project. The cellulose 
chain has actually a somewhat flattened shape 
and therewith two [2] hydrophobe and two [2] 
hydrophile sides.{ With no other hexose than 
8-glucose can such a relation be obtained. 

It is apparent that the configuration suggested 
by these authors would involve somewhat dif- 
ferent conditions for the adsorption of water 
vapor, and for its effect upon electrical anisot- 
ropy, from those deducible from the earlier 
Meyer-Mark chain macromolecule. However, 
evidently in one important respect they are 
alike, viz., in predicating a duplex or twofold 
array of adsorption ‘‘sites’’ for water to form 
continuous H-bond chains. The temperature 
coefficient of adsorption, or, more precisely, the 
differential adsorption heat, would seem likely 
to be of different magnitude and trend for the 
two cases, and may help to decide between them. 
The newer model, especially if the preferred 
formation of H bridges between glucosan —OH 
groups and bridge oxygens be accepted, would 
appear to require longer chains of oriented ad- 
sorbed water molecules than the earlier one. 

We may note in this connection that S. Baxter® 
has suggested that the electrical conduction of 
textile fibers, e.g., wool, involves a mechanism 
other than electrolytic. He suggests an electronic 
conduction, as in a semiconductor. It is assumed 
that ‘‘the dry wool fiber is an almost perfect 
insulator, and that adsorbed water molecules 
cause ‘impurity centers’ giving rise to electronic 
conduction.’ We are not convinced, for reasons 
which will be discussed elsewhere,** that much 
of the direct (macro) conductance is not elec- 
trolytic, but agree in part with Baxter for an 
inner, residual electronic portion. The actual 
mechanism of this we consider to be by way of a 
hydrogen-bond resonance process, as indicated 
here. 

A study of the temperature coefficients is con- 
templated, as well as of other factors, in the 
relation of electrical anisotropy to electrical 
conductance. 


t Our italics. 
*S. Baxter, Trans. Faraday Soc. 39, 207 (1943). 
** In connection with the conductance of gelatin. 
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Photography of Crystal Structures 


Maurice L. HuGGIns 
Eastman Kodak Company, Rochester, New York 
November 1, 1944 


S shown by the writer! several years ago, Sir Lawrence 
Bragg’s? method of photographic summation of 
Fourier series by addition of the proper patterns of light 
and dark bands can be improved in accuracy and speed by 
the use of a previously prepared set of masks. 





Fic. 1. Photograph, looking along one of the crystallographic axes, of 
the body-centered cubic unit of a-iron, from x-ray data by Armstrong. 


With the assistance of the Physics Department of these 
Laboratories, an improved set of 316 masks for this purpose 
has recently been made. Using these, electron density 
photographs of crystal structures and Patterson summa- 
tions are easily and rapidly made from the appropriate 
x-ray data. Examples are shown in Figs. 1-3. The magnifi- 
cation is of the order of 100,000,000 in each case. 

Most of the background density, between the atoms, 
in Fig. 33 is not real, but is a result of inaccuracies in the 
x-ray data and in the masks and incompleteness of the 
former. The relative positions and approximate relative 
densities of the atomic peaks are correctly shown, however. 
In Fig. 1,4 the background was removed by a short treat- 
ment with Farmer’s reducer solution. In Fig. 25 the back- 
ground has been practically all removed by reprinting. 

The shape of the molecule in Fig. 2 is somewhat dis- 
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Fic. 2. Photograph of the hexamethylbenzene molecule, using x-ray 


data from the crystal, by Brockway and Robertson. 


torted, both because the plane of the molecule is not 
parallel to the plane of the projection and because the 
projection, as obtained, is a square, whereas the true pro- 
jection of the unit cell is not. The latter cause of distortion 
can readily be removed, if desired, by adding another 
optical step or in other ways. For purposes of structure 
analysis, however, this is not necessary. 





Fic. 3. cogs oy ny of a B-resorcinol (meta-dihydroxybenzene) crystal, 
using x-ray data by Robertson and Ubbelohde. Two complete mole- 
cules, with parts of others, are shown. 


We hope soon to be able to furnish duplicate copies of 
the improved set of masks, at a nominal cost, to others 
engaged in crystal structure analysis. 


1M. L. Huggins, J. Am. Chem. Soc. 63, 66 (1941). 
2W. L. Bragg, Zeits. f. Krist. A70, 475 (1929); The gang A State 
(The Macmillan Company, London and New York, 1934), p. 
(1938), M. Robertson and A. R. Ubbelohde, Proc. Roy. Soc. ‘A167, 122 
4A, H. Armstrong, Phys. Rev. 34, 931 (1929). 
5L.O. Brockway and J. M. Robertson, J. Chem. Soc., p. 1324 (1939). 
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Errata : On the Statistical Mechanics of Liquids, 
and the Gas of Hard Elastic Spheres 
[J. Chem. Phys. 12, 1 (1944) ] 
O. K. RICE 
University of North Carolina, Chapel Hill, North Carolina 
November 2, 1944 

AGE 4: Third and fourth lines above Eq. (3) and else- 
where on the page. The ‘‘dodecahedron”’ is actually a 
figure with twelve vertices but fourteen sides (six squares 

and eight triangles). The volumes are given correctly. 

Page 11: Eq. (24). Second term in bracket, 
1 o 1 o 

, should be - , 
2 r,3¢8 


2 73a 








Page 16: Left-hand side of displayed unnumbered equa- 
tion slrould be 
[e/87-/@ (av) (242/24), 


Page 18: In line following first equation in Appendix, 
‘free space’’ should read “portion of the total length.” 

Note: It may be remarked that the term “cubic close- 
packed” and derivatives thereof are used in two senses. 
Sometimes they refer to a cubic close-packed array of hard 
spheres in mutual contact, which thus occupy the smallest 
possible volume. In other cases they refer rather to the 
type of lattice, without demanding that the spheres b2 in 
actual contact, the possibility of displacements about the 
average positions which form the lattice points being 
visualized. The actual usage in any given instance should 
be clear from the context. 





Erratum: The Chain Photolysis of Acetaldehyde 
in Intermittent Light 
[J. Chem. Phys. 10, 445 (1942) ] 


W. L. HADEN, JR., AND O. K. RICE 
University of North Carolina, Chapel Hill, North Carolina 
November 2, 1944 


N the abstract the energies of activation for reactions 
(2) and (3) are interchanged. They are given correctly 
in Table IIT. 





Bond Moments of Higher Valence States 


R. SAMUEL 


The Illinois Institute of Technology, Chicago, Illinois, and the 
Hebrew Institute of Technology, Haifa, Palestine 


November 10, 1944 


- a recent article! additive bond moments have been 
calculated for non-aromatic first-order molecules formed 
by a central atom in its maximal state of valency. This has 
been done on the basis of classical formulae. One point 
apparently has not been made sufficiently clear, as it has 
subsequently been raised in a number of discussions. This 
is, whether the new values in reality are not those of the 
semipolar double bonds. In other words, whether, for ex- 
ample, the value of 4.4D ascribed to the bond between 
pentacovalent nitrogen and oxygen in the nitro group, 
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does not constitute in reality the value for the semipolar 
double bond between the same two atoms. In this letter it 
will be shown that this is not so. 

On the basis of the resonance interpretation of the semi- 
polar double bond, the nitro group is a hybrid of the two 
forms 





— 
O=N—O (1a) 
| 
and 
- + 
O—N=O. (1b) 





| 


On account of the rigorous energetical degeneracy each 
form contributes 50 percent. Indeed, if no other similar 
molecule existed, the value of 4.4D could easily be asso- 
ciated with this mixture. 

But in the azoxy compounds the two corresponding forms 
are 


—N=N—O (2a) 
| 
and 
- + 
—N—N=0O. (2b) 


| 
| 


These constitute two different terms of the molecule, and 
according to their energetical difference (2a) will contri- 
bute x: percent and (2b) y: percent. To obtain the value 
of 4.4D for the mixture, the energy difference must be 
such that x; and y; just have the required values. This 
would be an entirely fortuitous coincidence, as this energy 
of excitation depends on both parts of the molecule equally 
and not on the nitrogen-oxygen part alone. 

Again the N-ethers of oximes would be mixtures of x2 
percent of 


>C=N—O (3a) 
and ye percent of 


>C—N=0, (3b) 
| 


and if fortuitously the energy difference between these 
two terms is such that x2 and y2 obtain the required ratio 
the value of 4.4D may result for the moment of this bond. 
The same would be true of nitrous oxide, which would 
be a mixture of x3 percent of 





N=N—O (4a) 
and ys; percent of 
- + 
N=N=0. (4b) 


This is only half of the problem. A very great number 
of molecules are interconnected by having one or the other 
bond in common. Thus, the same bond between C and N 
of the N-ethers occurs again with the same moment in the 
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diazo compounds which would be mixtures of xs percent 
and ys percent of two forms. The ratio of x, and y, must 
be such as to obtain the same value of 2.6 D for the 
bond between C and N as in the N-ethers. But for the 
latter the percentage ratio of x2 and y2 has already been 
fixed according to the requirements of the bond between 
N and O. The difference between the two energy levels 
of the N-ethers and hence the ratio of x2 and y2 must be 
such as to obtain for this mixture of two forms the same 
moment of the nitrogen-oxygen bond as in the nitro group, 
of the nitrogen-carbon (double) bond the same as in the 
diazo compounds, and of the nitrogen-carbon (single) bond 
the same as in the nitrates. The argument must not be 
elaborated. But it is very difficult to see how the excita- 
tion energies of the molecules can be consistent with so 
many different and simultaneous conditions. 

To sum up: The problem which confronts us cannot be 
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solved by means of philology but is a genuine problem in 
the realm of physics. On the basis of the classical formulae 
we expect in first approximation constant and additive 
moments of the individual bonds, and we find them. On 
the basis of the resonance interpretation of the theory of 
the semipolar double bond we really cannot expect addi- 
tive and constant bond moments, but we find them 
nevertheless. Hence we have to explain each of them 
by a whole series of fortuitous coincidences, which in the 
framework of this theory remain unexplained freaks of 
nature. As has been pointed out, such a state of affairs 
cannot be directly disproved by a physico-chemical experi- 
ment, because the percentage ratio is a variable factor in 
the latter theory. But it would appear that the classical 
concepts are superior by virtue of the principle of economy 
of hypotheses. 


1R. Samuel, J. Chem. Phys. 12, 380 (1944). 
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